Dark energy perturbations and cosmic coincidence 
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While there is plentiful evidence in all fronts of experimental cosmology for the existence of a 
non-vanishing dark energy (DE) density fjn in the Universe, we are still far away from having a 
fundamental understanding of its ultimate nature and of its current value, not even of the puzzling 
fact that pd is so close to the matter energy density pM at the present time (i.e. the so-called 
"cosmic coincidence" problem). The resolution of some of these cosmic conundrums suggests that 
the DE must have some (mild) dynamical behavior at the present time. In this paper, we examine 
some general properties of the simultaneous set of matter and DE perturbations {5pM , Spn) for 
a multicomponent DE fluid. Next we put these properties to the test within the context of a 
non-trivial model of dynamical DE (the AXCDM model) which has been previously studied in the 
literature. By requiring that the coupled system of perturbation equations for SpM and Spn has a 
smooth solution throughout the entire cosmological evolution, that the matter power spectrum is 
consistent with the data on structure formation and that the "coincidence ratio" r = po/pM stays 
bounded and not unnaturally high, we are able to determine a well-deflned region of the parameter 
space where the model can solve the cosmic coincidence problem in full compatibility with all known 
cosmological data. 

PACS numbers: 95.36.-l-x, 04.62.-(-v, ll.lO.Hi 



I. INTRODUCTION 



Undoubtedly the most prominent accomplishment of 
modern cosmology to date has been to provide strong 
indirect support for the existence of both dark matter 
(DM) and dark energy (DE) from independent data sets 
derived from the observation of distant supernovae [l[ , 
the anisotropics of the CMB 2] , the lensing effects on the 
propagation of hght through weak gravitational fields 3] , 
and the inventory of cosmic matter from the large scale 
structures (LSS) of the Universe 0, But, in spite 
of these outstanding achievements, modern cosmology 
still fails to understand the ultimate physical nature 
of the components that build up the mysterious dark 
side of the Universe, most conspicuously the DE compo- 
nent of which the first significant experimental evidence 
was reported 10 years ago from supernovae observations. 
The current estimates of the DE energy density yield 
~ (2.4 X 10-3 eV)^ and it is believed that it con- 
stitutes roughly 70% of the total energy density budget 
for an essentially flat Universe. The big question now is: 
what is it from the point of view of fundamental physics? 
One possibility is that it is the ground state energy den- 
sity associated to the quantum field theory (QFT) vac- 
uum and, in this case, it is traditional to associate pn to 
PA = A/SttG, where A is the cosmological constant (CC) 
term in Einstein's equations. The problem, however, is 
that the typical value of the (rcnormalized) vacuum en- 
ergy in all known realistic QFT's is much bigger than 



the experimental value. For example, the energy density 
associated to the Higgs potential of the Standard Model 
(SM) of electroweak interactions is more than fifty orders 
of magnitude larger than the measured value of pu. 

Another generic proposal (with many ramifications) is 
the possibility that the DE stands for the current value of 
the energy density of some slowly evolving, homogeneous 
and isotropic scalar field (or collection of them) . Scalar 
fields appeared first as dynamical adjustment mecha- 
nisms for the CC [1,0| and later gave rise to the notion of 
quintessence Q . While this idea has its own merits (spe- 
cially concerning the dynamical character that confers to 
the DE) it has also its own drawbacks. The most obvious 
one (often completely ignored) is that the vacuum energy 
of the SM is still there and, therefore, the quintessence 
field just adds up more trouble to the whole fine-tuning 
CC problem I, [l|! 

Next-to-leading is the "cosmological coincidence pro- 
blem" , or the problem of understanding why the 
presently measured value of the DE is so close to the mat- 
ter density. One expects that this problem can be allevi- 
ated by assuming that is actually a dynamical quan- 
tity. While quintessence is the traditionally explored op- 
tion, in this paper we entertain the possibility that such 
dynamics could be the result of the so-called cosmologi- 
cal "constants" (like A, G,...) being actually variable. It 
has been proven in [ll| that this possibility can perfectly 
mimic quintessence. It means that we stay with the A 
parameter and make it "running" , for example through 
quantum effects [H [H [ll ^ However, in [17] it was 
shown that, in order to have an impact on the coincidence 
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problem, the total DE in this context should be con- 
ceived as a composite fluid made out of a running A and 
another entity X, with some effective equation of state 
(EOS) parameter wx, such that the total DE density and 
pressure read po = Pa + px and pd = —pA+i^xPx, res- 
pectively. We call this system the AXCDM model [13] ■ 
Let us emphasize that X (called "the cosmon") is not 
necessarily a fundamental entity; in particular, it need 
not be an elementary scalar field. As remarked in p7j . 
X could represent the effective behavior of higher order 
terms in the effective action (including non-local ones). 
This is conceivable, since the Bianchi identity enforces 
a relation between all dynamical components that enter 
the effective structure of the energy-momentum tensor 
in Einstein's equations, in particular between the evolv- 
ing A and other terms that could emerge after we embed 
General Relativity in a more general framework [l^, [l^ . 
Therefore, at this level, we do not impose a microscopic 
description for X and in this way the treatment becomes 
more general ^ . The only condition defining X is the DE 
conservation law, namely we assume that pd — PA Px 
is the covariantly self-conserved total DE density. 

In this paper, we analyze the combined dynamics of DE 
and matter density perturbations for such conserved DE 
density po- The present stu dy g oes beyond the approxi- 
mate treatment presented in [23] , where we neglected the 
DE perturbations and estimated the matter perturba- 
tions of the AXCDM model using an effective (variable) 
EOS We for the composite fluid {pA,px)- The main re- 
sult was that a sizeable portion of parameter space was 
still compatible with a possible solution of the cosmic 
coincidence problem. The "effective approach" that we 
employed in was based on three essential ingredi- 
ents: i) the use of the effective EOS representation of 
cosmologies with variable cosmological parameters [ll[ ; 
ii) the calculation of the growth of matter density fluc- 
tuations using the effective EOS of the DE [23|; and iii) 
the application of the, so-called, "F-test" to compare the 
model with the LSS data, i.e. the condition that the lin- 
ear bias parameter, iP'^z) = Pgg/Pmm does not deviate 
from the ACDM value by more than 10% at z = 0, where 
Pmm oc {5pM / PmY is the matter power spectrum and 
Pan is the galaxy fluctuation power spectrum [1, Q - see 
[23I . [2BI for details. This three-step methodology turned 
out to be an efficient streamlined strategy to further 
constrain the region of the original parameter space . 
However, there remained to perform a full fledged analy- 
sis of the system of cosmological perturbations in which 
the DE and matter fluctuations are coupled in a dynam- 
ical way. This kind of analysis is presented here. 

The structure of the paper is as follows. In the next 
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section, we outline the meaning of the cosmic coincidence 
problem within the general setting of the cosmological 
constant problem. In section Hill the basic equations for 
cosmological perturbations of a multicomponent fluid in 
the linear regime are introduced. In section IIVI we de- 
scribe the general framework for addressing cosmologi- 
cal perturbations of a composite DE fluid with an effec- 
tive equation of state (EOS). In section |Vl we describe 
some generic features of the cosmological perturbations 
for the dark energy component. The particular setup of 
the AXCDM model is focused in section IVIl In sections 
rVlIland lVini we put the AXCDM model to the stringent 
test of cosmological perturbations and show that the cor- 
responding region of parameter space becomes further re- 
duced. Most important, in this region the model is com- 
patible with all known observational data and, therefore, 
the AXCDM proposal can be flnally presented as a ro- 
bust candidate model for solving the cosmic coincidence 
problem. In section HXl we offer a deeper insight into the 
correlation of matter and DE perturbations. In the last 
section, we present the final discussion and deliver our 
conclusions. 



II. THE COINCIDENCE PROBLEM AS A PART 
OF THE BIG CC PROBLEM 

The cosmic coincidence problem is a riddle, wrapped 
in the polyhedric mystery of the Cosmological Constant 
Problem [l,[l3], which has many faces. Indeed, we should 
clearly distinguish between the two main aspects which 
are hidden in the cosmological constant (CC) problem. In 
the first place, we have the "old CC problem" (the ugliest 
face of the CC conundrum!), i.e. the formidable task of 
trying to explain the relatively small (for Particle Physics 
standards) measured value of pA or, more generally, of 
the DE density roughly ~ 10"^^ GeV^, after the 
many phase transitions that our Universe has undergone 
since the very early times, in particular the electroweak 
Higgs phase transition associated to the Standard Model 
of Particle Physics, whose natural value is in the ballpark 
oi Pew ^ Gp^ 10^ GeV^ {Gp being Fermi's constant). 
The discrepancy Pew / Pd^ : which amounts to some 56 
orders of magnitude, is the biggest enigma of fundamen- 
tal physics ever! Apart from the induced CC contribu- 
tion from phase transitions, we have the pure vacuum-to- 
vacuum quantum effects. Since the (renormalized) zero 
point energy of a free particle of mass m contributes 
^ to the vacuum energy density [l^, [l^ , it turns out 
that even a free electron contributes an amount more 
than thirty orders of magnitude larger than the afore- 
mentioned experimental value of pu. Only light neutri- 
nos nil, ~ 10"'^ eV, or scalar particles of similar mass. 
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could contribute just the right amount, namely if these 
particles would be the sole active degrees of freedom in 
our present cold Universe (see [l^). 

On the other hand, the cosmic coincidence problem f2^ 
is that second ("minor") aspect of the CC problem ad- 
dressing the specific question: "why just now?" , i.e. why 
do we find ourselves in an epoch t — to where the DE 
density is similar to the matter density, poito) — Pm(^o)? 
In view of the rapidly decreasing value of Pm{o) ^ l/a"^, 
where a = a{t) is the scale factor, it is quite puzzling 
to observe that its current value is precisely of the same 
order of magnitude as the vacuum energy or, in general, 
the dark energy density po . It is convenient to define the 
"cosmic coincidence ratio" 



r(a) 



PAiia) ^M{a) ' 



(1) 



where (r2D(a), i7M(a)) are the corresponding densities 
normalized with respect to the current critical density 
p° EE 3H^/8nG. For ~ 0.3 and ~ 0.7, we have 
rg ~ 2.3, which is of 0(1). However, in the standard 
cosmological ACDM model, where Jl^i is constant and 
^Mia oo) 0, the ratio r grows unboundedly with 
the expansion of the Universe. So the fact that tq — 0{1) 
is regarded as a puzzle because it suggests that t — to is 
a very special epoch of our Universe. One could also con- 
sider the inverse ratio = Pm{<^) I Pd{o)-, which goes to 
zero with the expansion. The coincidence problem can 
be equivalently formulated either by asking why r is not 
very large now or why is not very small. Solving the 
coincidence problem would be to find either 1) a concrete 
explanation for r and being of order one at present 
within the standard cosmological model, or 2) a modified 
cosmological model (compatible with all known cosmolo- 
gical data) insuring that these ratios do not undergo a 
substantial change, say by more than one order of mag- 
nitude or so, for a very long period of the cosmic history 
that includes our time. 

In a very simplified way, let us summarize some of the 
possible avenues that have been entertained to cope with 
the coincidence puzzle: 

• Quintessence and the like [1, [13, [H, [H, HO] ■ One 

postulates the existence of a set of cosmological 
scalar fields (pi essentially unrelated to the rest of 
the particle physics world. The DE produced by 
these fields has an effective EOS parameter lod ^ 
— 1 which causes to decrease always with the 
expansion (i.e. dpi:) /da < 0), but at a pace slower 
(on average) than that of the background matter. 
Thus, it finally catches up with it and pjj emerges 
to surface, i.e. the condition po > Pm eventually 
holds (presumably near our time). In this frame- 
work, there is the possibility of self-adjusting and 
tracker solutions [1, [13] , where the DE keeps track 
of the matter behavior and ultimately dominates 
the Universe. It requires to take some special forms 
of the potential, and in some cases the Lagrangian 



involves non-canonical kinetic terms. For example, 
in the simple case of a single scalar field and the 
exponential potential V{(j)) ~ exp {—X(j)/Mp) one 
finds that the coincidence ratio becomes fixed at 
the value 



3(1 



A2 -3(l+a;„) 



(2) 



where LOm is the EOS of the background matter 
(i.e. or 1/3, depending on whether cold or rel- 
ativistic matter dominates, respectively). So, at 
the present time, r — 3/(A^ — 3), and by appro- 
priate choice of A one can match the current ex- 
perimental value. But of course the choice of the 
potential was rather peculiar and the field (j) it- 
self is completely ad hoc. Moreover, it has a mass 
= y/V^ H ^ 10-33 (as it follows 
from a self-consistent solution of Einstein's equa- 
tions); such mass scale is 30 orders of magnitude 
below the mass scale associated to the DE, which 
is p}/^ ~ 10"'^ eV . In this sense, it looks a bit 
unnatural to aim at solving the CC problem by in- 
troducing a field whose extremely tiny mass creates 
another cosmological puzzle. 

On the other hand, within the context of interac- 
tive quintessence models [28i] (whose main leitmo- 
tif is precisely trying to cure the coincidence pro- 
blem), the coupling of 4)i and the matter compo- 
nents makes allowance for energy exchange between 
the two kinds of fields, and as a result the ratio 
([1]) can be constant or slowly variable, whereas in 
other implementations one can achieve an oscilla- 
tory tracking behavior of r, although the construc- 
tion is essentially ad hoc [29.] . Another generaliza- 
tion leads to k-essence models (characterized 
by non-canonical kinetic terms), where fine tuning 
problems in the tracking can be disposed of but the 
dominant background component can be tracked 
only up to matter-radiation equality and is lost 
immediately afterwards (as the DE is immediately 
prompted into a CC-like behavior). In one way or 
another, however, all variants of quintessence suf- 
fer from several drawbacks, and in particular from 
the following generic one: they assume (somehow 
implicitly) that the remaining fields of the parti- 
cle physics spectrum (i.e. those which were already 
there from the very beginning) have nothing to do 
with the CC problem. As a result of such a bold as- 
sumption, the (likely real) vacuum problem of the 
conventional fields in QFT is merely traded for the 
(likely fictitious) vacuum problem of quintessence, 
which is no less acute because no real explanation is 
provided for the smallness of the current value 
versus (where m is any typical mass scale in 
Particle Physics). Hence we are back to the same 
kind of CC problem we started with. 

• Phantom energy [3l[ . It is motivated by the fact 
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that, observationally speaking, the effective EOS 
of the DE cannot be excluded to satisfy ujjj < —I 
near our present time. As indicated above, many 
quintessence-hke models in reality are hybrid con- 
structions containing a mixture of fields with a 
phantom component. The reason is that one wants 
to give allowance for a "CC-crossing" uiu = —1 
near our time. While phantom energy shares with 
quintessence the use of scalar fields here the DE 
produced by these fields is always increasing with 
the expansion, dpu/da > 0, even after the relation 
Pd > Pm is fulfilled. The consequence of this ever 
growing behavior of the dark component is that one 
ends up with a superaccelerated expansion of the 
Universe that triggers an eventual disruption of all 
forms of matter (the so-called "Big Rip"). When 
computing the fraction of the lifetime of the Uni- 
verse where the ratio ([T]) stays within given bounds 
before the "doomsday", one finds that it can be 
sizeable. 

• Non-local theories. There is some renewed interest 
in this kind of theories in which the emphasis is 
placed on the existence of possible non-local struc- 
tures in the effective action ^^l- It has recently 
been emphasized in [l6l| that the dynamical evolu- 
tion of the vacuum energy should come from a re- 
summation of terms in the effective action leading 
to non-local contributions of the form RT{QqR), 
for some unknown function J- of dimension 2, where 
Qq is the massless Green's function (C/o ~ l/D)- 
The canonical possibility would he T — APGqR, 
where M is a parameter with dimension of mass. 
This situation leads to an effective evolution of the 
CC of the form ApA ~ A/^ during the matter 
dominated epoch, whereas in the radiation era the 
effective CC would approximately be zero (because 
i? ~ ~ for relativistic matter, see (fT2t - (fT3)) 
below). As a result, the coincidence puzzle could 
somehow be understood from the fact that the CC 
may start to be preponderant at some point once 
the onset of the matter dominated epoch is left be- 
hind. 

• Of course many other ideas have been explored. For 
instance, one may introduce special fluids with very 
peculiar EOS, such as the Chaplygin gas [33|, which 
behaves as pressureless matter at early times {ujd — 
0) and as vacuum energy at present {ljd — —1). Al- 
though there is some connection with braneworld 
cosmology, this proposal suffers from the same pro- 
blem as quintessence in that it supersedes the vac- 
uum state of traditional flclds by the new vacuum of 
that peculiar fluid. Finally, let us mention the An- 
thropic models, which fall in a quite different cate- 
gory, in the sense that one does not look for a so- 
lution of the coincidence problem exclusively from 
first principles of QFT or string theory, but rather 
through the interplay of the "human factor" . Basi- 



cally, one ties the value of the ratio ([T]) to the time 
when the conditions arise for the development of 
intelligent life in the Universe, in particular of cos- 
mologists making observations of the cosmos. This 
variant has also a long story, but we shall refrain 
from entering the details, see e.g. [ol. [33|. 



III. COSMOLOGICAL PERTURBATIONS FOR 
A MULTICOMPONENT FLUID IN THE LINEAR 
REGIME 

In the remaining of the paper we concentrate on study- 
ing some general properties of the cosmological pertur- 
bations, both of matter and DE, and the implications 
they may have on the coincidence problem within mo- 
dels characterized by having running vacuum energy and 
other DE components. According to cosmological per- 
turbation theory, all energy density components, includ- 
ing the dark energy, should fluctuate and contribute to 
the growth of the large-scale cosmological structures. In 
this section, we discuss the general framework of linear 
density perturbations in models composed of a multicom- 
ponent DE fluid besides the canonical matter. 

In the following we use the standard metric perturba- 
tion approach [35| and consider simultaneous density and 
pressure perturbations 

Pn ^ Pn + SpN , Pn ^ Pn + SpN , (3) 

for all the components {N — 1,2,...) of the fluid, in- 
cluding matter and all possible contributions from the 
multicomponent DE part. At the same time, we have 
metric and velocity perturbations for each component: 

g!^. ^ 9^. = g!^. + , u^-^ui:, + su^ , (4) 

where g^^^ is the background metric. The 4-vector ve- 
locity [/^ in the comoving coordinates has the follow- 
ing components and perturbations (t/^,t/^) = (1,0) 
6U^ — (0,w]y), where v]^ is the three-velocity of the 
iVth component of the fluid in the chosen coordinate sys- 
tem. As a background space-time, we assume the homo- 
geneous and isotropic Friedmann-Lemaitre-Robertson- 
Walker (FLRW) metric with flat space section, hence 
.9^1^ = diagjl, —a^{t) 6ij}, where a is the scale factor. 

In order to derive the set of perturbed equations, let 
us first introduce Einstein's equations: 

Rf^iy ~ ^Rgfj.iy ^ 8ttGT^^ , (5) 

where G is the Newton constant and T^^ is the total 
energy-momentum tensor of matter and dark energy. 
Both the background and perturbed metric are assumed 
to satisfy these equations. The total energy- momentum 
tensor of the system is assumed to be the sum of the 
perfect fluid form for each component: 

Tf.. = E ^M- = E [-P^ + {pN + PN ) U^] . 
N N 

(6) 
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The components of are the following: 

T^-^PAT^PT, (7) 

N 

T} = -Y.P^S]^-prS], (8) 

N 

where pT and pr are the total energy density and pres- 
sure, respectively. 

Perturbations on the metric and on the energy- 
momentum tensor are uniquely defined for a given per- 
turbed space-time, provided we make a gauge choice. 
The latter means that we choose a specific coordinate 
system; in this way, four out of the 10 components 
Sgfj^u = hfj^i, of the metric perturbation can be fixed 
at will. Here we have adopted the synchronous gauge, 
widely used in the literature, in which the four preas- 
signed values of the metric perturbations are /loi = 
and /loo = 0. Setting h^^i, = — a^x^^, in this gauge the 
perturbed, spatially flat, FLRW metric takes on the form 

ds^ ~ Qiiv dx^dx^ = dt^ — a^{t) (Sij + Xij) dx^dx-' 
= 0^(77) [dT^^ - -f Xij) dx'dx'] , (9) 

where in the last equality we have expressed the result 
also in terms of the conformal time 77, defined through 
drj = dt/a. We may compare ^ with the most general 
perturbation of the spatially flat FLRW metric, 

ds^ = a^{r])[{l + 2i}})dif -ujidi^dx' 

- {{l-2d^)5,J+x^3)dx'dx^], (10) 

consisting of the 10 degrees of freedom associated to the 
two scalar functions ip^ 4>, the three components of the 
vector function uti (i — 1,2,3), and the five components 
of the traceless Xij ■ Clearly, the synchronous gauge ([9]) 
is obtained by setting = 0, cui — and absorbing the 
function (p into the trace of Xij ■ In this way, Xij i'^ ® 
contributes six degrees of freedom. As we will see in a 
moment, in practice only the nonvanishing trace of the 
metric disturbance will be necessary to perform the anal- 
ysis of cosmic perturbations in this gauge. To within first 
order of perturbation theory, such a trace is given by 

/i^g^^'V-s'-' /i^. --^=X.», (11) 

where g^^ is the inverse of gij = — a^(i) {Sij+Xij)i and it is 
understood that the repeated Latin indices are summed 
over 1, 2, 3. 

For the physical interpretation, notice that the syn- 
chronous gauge is associated to a coordinate system in 
which the cosmic time coordinate is comoving with the 
fluid particles (500 = 1, i-e. V = 0), which is the rea- 
son for its name and also explains why this gauge does 
not have an obvious Newtonian limit. In fact, this gauge 
choice is generally appropriate for the study of fluctua- 
tions whose wavelength is larger than the Hubble radius 



(A ^ dH = H~^). Actually, any mode satisfies this con- 
dition at sufficiently early epochs, and in this regime the 
effects of the space-time curvature are unavoidable. 

Next we wish to compute the perturbations in the syn- 
chronous gauge. To start with, it is convenient to rewrite 
Einstein's equations ^ as follows, 

Rfj.v — SttG S 7 S = T^i, — g^i, T , (1^) 

where T — Tj^ \s the trace of ([S]), hence 

T^ = II(Pw-3pw) . (13) 

N 

For the calculation of the perturbations, we can use any 
of the components of Einstein's equations. However, 
since we are going to use the conservation law V^T^ = 
to derive additional fluctuation equations, it is conve- 
nient to perturb the (OO)-component of (fT2|) because the 
other components are well-known not to be independent 
of the conservation law. Thus, using ©, (fT^ and (fT5|) 
we obtain 

•^oo = ^00 - ^ 7" = ^ ^ {PN + 3pAr) 

N 

N 

On the other hand, a straightforward calculation shows 
that the perturbation of the (OO)-component of the Ricci 
tensor can be written as follows: 

1 d'^h dh I ° 
where we have used ((TT]) and defined the "hat variable" 

The overhead circle ( ° ) indicates partial differentiation 

o 

with respect to the cosmic time (i.e. /= df /dt, for any 
/), in order to distinguish it from other differentiations 

to be used later on. Therefore, H = a /a is the ordinary 
expansion rate in the cosmic time t. 

Since the fluctuations ^6*00 and SRqq from (|14p and 
(fT5|) are constrained by (fT2|) . we flnd 

o 

h+2Hh^8'KGj2i^PN + ^SpN). (17) 

N 

If we substitute (|16p in the previous expression, a second 
order differential equation in the original variable (jlip en- 
sues. In terms of the conformal time 77, it can be written 
as follows: 

h + nh = -SttGo^ ^{SpN + 3Spn) , (18) 

N 
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where the dot ( ' ) indicates differentiation with respect to 
r] (i.e. / = df / drj) and T-l = a/a = Ha is the expansion 
rate in the conformal time. 

The Friedmann equation can be written in terms of the 
normaUzed densities as 



-PT 



Hi ^ r2jv 



(19) 



N 



where Hq is the present value of the Hubble parameter 
and rjjv = Pn / p'i are the normaUzed densities with re- 
spect to the current critical density = ZH^ / 8ttG. 

The subsequent step is to perform perturbations on the 
conservation law V^T^'' = 0, as this will provide the ad- 
ditional independent equations. Using ([S]), the previous 
law reads explicitly as follows: 

J2 {^m(pw + Pn) U^U^ + {pN + Pn) 

N 

N 

(20) 

For any four-velocity vector, we have U'^ Uj^ — 1 and, 
therefore, we have the orthogonality relation U^V ^JJ^ — 
0. In this way, by contracting Eq. (PO)) with Uj^ we find 
the simpler result 



^ [C/j^ PN + {pN + PJV)V^[/^] = . 



(21) 



N 



Let us emphasize that the sum in this equation need 
not run necessarily over all the terms of the cosmic fluid. 
It may hold for particular subsets of fluid components 
that are overall self-conserved. In particular, it could 
even hold for each component, if they would be individ- 
ually conserved. In our case, it applies to the specific 
matter component and also, collectively, to the multi- 
component DE part. It is straightforward to check that, 
in the FLRW metric, we have 



V^U^ = 3H (TV = 1,2,...). 



(22) 



Using this relation, it is immediate to see that, in the 
co-moving frame, Eq. l|2T|) boils down to 



pAT +3H{pn+Pn) 

N 



= 0. 



(23) 



Moreover, perturbing (|22p in the synchronous gauge, we 
find (using ST'^^ = —h/2 for the perturbed Christoffel 
symbol involved in the covariant derivative) the useful 
result 



On 



(24) 



where we have introduced the notation 6 
Vr{5Ul,) (with5[/o 



jv = V^(5(7^^) = 
0) for the covariant derivative of 



the perturbed three-velocity 5Uli 
these formulas, the perturbed Eq. 
to 



]y. Equipped with 
immediately leads 



E 

N 



dpN +ipN + Pn) \(^n - ] + 3H{dpN + Spn) 



^0. 



(25) 

The previous result could have equivalently been ob- 
tained by setting = in (PO)) and perturbing the cor- 
responding equation. An independent relation can be 
obtained by setting v — i in (|20p within the co-moving 
frame and carrying out the perturbation. Using the rel- 



evant Christoffel symbol Fq^- 
tain, after some calculations. 



HSij and Eq.dm we ob- 



N 



E 1 (/^^ + ^^'n + \°Pn +Pn +5H{pN + pn) 



X SU 



N 



Y^y'SpN. 



(26) 



N 



The final step is obtained by computing the divergence 
Vi on both sides of this equation. To within first 
order of perturbation theory, we obtain "S/iS/^SpN — 
g^''\/i\/kSpN = —{l/a?)y^5pN for the action of Vi on 
the r.h.s. of (|26|) . whereas on its l.h.s. we can use 
the variable On previously defined. Moreover, we have 
° ,■ ° . 

Vi 5U'^ N ^ 9n owing to the commutativity of coordinate 
differentiation and perturbation operations. In this way, 
the final outcome reads 

YI^P^ + On + [°Pn +Pn +5H{pN + pn) 



N 



oa='4j:^p 



N, 



(27) 



N 



where it is furthermore understood that we have used the 
Fourier decomposition for all the perturbation variables: 



(5/(t,x) 



(27r)^ 



6f{t,k)e'^--. 



(28) 



where Sf = {h,5pN,5pN,()N)- In Fourier space, the per- 
turbation variables are denoted with the same notation, 
but they are the Fourier transforms of the original ones, 
so their arguments are t and k because the space variable 
X has been traded for the wave number k= |k|. The lat- 
ter will be measured in units of /iMpc~^, where h ~ 0.7 
is the reduced Hubble constant - not to be confused with 
the trace of the synchronous perturbed metric, Eq. (jlip . 
In these units, the linear regime corresponds to length 
scales £ ^ fc^^ with wave numbers k < 0.2ft,Mpc~^, i.e. 
£ > bh~^ Mpc. Notice that, if desired, one can easily 
rewrite the above perturbation equations ([25]) and ([27|) 



in conformal time simply by using /= //a for any /. 

We have obtained three basic sets of perturbation 
equations p?|) . (^5)) and (P7)) for the four kinds of pertur- 
bation variables (/i(t, fc), 6pN(t,k), SpN{t,k), 9N{t,k)). 
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It is thus clear that the evolution of the cosmic pertur- 
bations can be completely specified only after we assume 
some relation of the pressure perturbation 5pN and the 
density perturbation 5pN for each fluid. If the perturba- 
tions are adiabatic, then that relation is simply 



= cI nSpn , 



(29) 



where is the adiabatic speed of sound for each fluid, 
defined as: 



^2 



Pn 

PN 



(30) 



where the dot differentiation here is with respect to what- 
ever definition of time. Notice that if the various compo- 
nents would have an equation of state (EOS) of the form 
Pn — wn Pn, with constant EOS parameter wn, then 
N ~ However, even in this case the mixture has 
a variable effective EOS parameter, as we will see in the 
next section. 

On the other hand, if the perturbations are non- 
adiabatic, there is an entropy contribution to the pressure 
perturbation [s^ : 



PnTn = 5pN - c^.n^Pn , 



(31) 



where Tn = ((5pAr)noii-adiabatic/PAr is the intrinsic en- 
tropy perturbation of the TVth component, representing 
the displacement between hypersurfaces of uniform pres- 
sure and uniform energy density [37| . For covariantly 
conserved components, a gauge-invariant relationship be- 
tween 5pN and 5pN for a general non- adiabatic stress is 
given by [Sj-H^: 



8pN = cl n5pn+ iHa^{pN+PN){c, 



,2 

s,N ^ 



where jy can be regarded as a rest frame speed of 
sound. We will refer to as the effective speed of 
sound, in the sense that we treat the cosmic fluid effec- 
tively as hydrodynamical matter. Since ((32l) is gauge 
invariant, the perturbed quantities in this expression 
can be computed, in particular, within the synchronous 
gauge. In this way, we can consistently substitute (|32p in 
the equations pT|) , (^5]) and (P7)) to eliminate the pertur- 
bation 5pN- This allows us, finally, to solve for the three 
basic sets of perturbations 



(/i(i, fc), SpN{t,k), dN{t,k)). 



(33) 



IV. 



PERTURBATIONS FOR A COMPOSITE DE 
FLUID WITH A VARIABLE EFFECTIVE 
EQUATION OF STATE 



In this section, we apply the linear matter and dark en- 
ergy density perturbations to a general class of models in 



which the DE fluid is a composite and covariantly self- 
conserved medium and matter is also canonically con- 
served. From Eq. ([^5]) . in the matter dominated epoch 
(pm = 0), the matter component p]\i satisfies 



Pnia) + -PM{a) = 0. 
a 



(34) 



Here we found convenient to trade the differentiation 
with respect to the cosmic time (°) for the differentiation 
with respect to the scale factor. The latter is denoted by 

o 

a prime (i.e. /' = df /da for any /; hence f— aH f). 
The scale factor is related to the cosmological redshift z 
by a{z) = 1/(1 + z), where we define a(0) = ao = 1 at 
the present time. It follows that the normalized matter 
density evolves as 



r^Af (a) 



where $7^^ is the corresponding current value, 
total DE is also globally conserved, from Eq. 
obtain 



p'd 



3 

+ -{l+We)pD ^0, 

a 



(35) 

Since the 
we also 



(36) 



where We is the effective equation of state (EOS) param- 
eter and pd = Pi + P2 + ■■■ is the total density of the 
multicomponent DE fluid. For a composite DE model, 
in which the DE is a mixture of fluids with individual 
EOS Pi — LOi Pi (i ~ 1, 2, n), the effective EOS param- 
eter is defined as 



Pd _ ^1 Pi +UJ2P2 + 
Pd P1+P2 + --- 



(37) 



The Hubble expansion in terms of the normalized den- 
sities, in the matter dominated period, follows from (|19|1 : 

= ^PT = nl {^\, a-' + Qnia)) , (38) 

where ilij is the normalized total DE density fl]j{a) = 
PD{a)/p°. 

The non-adiabatic perturbed pressure ([5^ for the total 
DE component can be written in terms of the effective 
EOS as 4 



SpD = C^SpD + 3i?a^(l + We)pD{c1 - <?a) 



fc2 



(39) 



where we have omitted for simplicity the subindex 'D' 
from the adiabatic and effective speeds of sound of the 
DE (i.e. = £1, = c^^; this convention will be 
used throughout the text). The units are taken to be the 



* In the particular case where the DE is a purely running A, one 
has to consider the perturbations 5pAi but since oja = — 1 it 
turns out that 6^3 (here 6a) plays no role and hence it is enough 
to consider the relation (5pA = ~^PKi see Ref. O for details. 
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light speed c = 1 and h = 1 such that the Planck scale 
is defined by Mp = G^i/^ ^ i 22 x 10^^ GeV. In these 



where in the last equation we have used pOj) to eliminate 
c^. Moreover, from ((45|) one can see that a negligible DE 



units, usually < < 1 for a general DE model. In this sound speed (c^ ~ 0) causes the velocity gradient to de- 



range, one can show that for constant EOS parameter 
there is small suppression on the DE fluctuation Spo as 
increases jSTj ] . Near-zero (but not vanishing) sound speed 
today is possible in models like k-essence, for example, 
in which the EOS parameter is positive until the matter- 
dominated triggers a change to negative pressure; in this 
kind of models, it is even possible to have > 1, regime 
for which t he g rowth of the DE density perturbations is 
suppressed |42 |. 

In a non-perfect fluid, spatial inhomogeneities in T^^ 
imply shear viscosity in the fluid. In this case, a possible 
contribution to shear through a "viscosity parameter" 



. In principle, 
the limit where 



c^ig should also be taken into account [43j 
is an arbitrary parameter. Nevertheless, 
(c^,c^is) (liO) corresponds exactly to a scalar field 
component with canonical kinetic term [4^ . 

Using the total DE conservation law ((36|) , we can write 
the total DE adiabatic sound speed ([30|) as 



Er 

p'd 



We 



3 {l + We) 



(40) 



The perturbed equations ([^S]) and ([7f|) for the (con- 
served) matter component (for which pM = SpM — 0) 
can be written as differential equations in the scale fac- 
tor: 



S'm = 



'M 




(41) 
(42) 



where 5m = 5pmI Pm is the relative matter fluctuation 
(density contrast). According to Ea. (|^ . the matter ve- 
locity gradient is decaying {6m oc a 
conventional initial condition = 
have 9m {a) — 0(Va). So, the perturbed matter set of 
coupled equations (HT1) and (|^^ yields the simple rela- 
tion 



^). Assuming the 
M{a = 1) = 0, we 



h = 2aHS'M ■ 



(43) 



Let us also define the relative fluctuation of the DE 
component, Sd = Spo/pD- Using the non- adiabatic per- 
turbed pressure ([55)1 and the DE conservation law 
we can write the perturbed equations and ([TT]) for 
the self-conserved DE fluid in the following way 



Si, 



aH 

3 

--{cl - We)S_ 

a 



9a'HHcl-cl) 



fc2 



D ■ 



(44) 
(45) 



cay (Od (X. a"^), as in the case of matter [Eq. (P^ ]. If we 
assume the conventional initial conditions = O'^ — 0, 
we have Om = do = (Va). In this case, the total DE 
fluid is comoving with the matter as long as the Uni- 
verse and perturbations evolve, which is a very partic- 
ular case. Actually, for this case, the k (scale) depen- 
dence disappears from the equations. On the other hand, 
from Eg. ([IS]) one can see that, if we neglect the DE per- 
turbations, SpD ~ 0, for a constant we obtain again 
^M = Od = and the scale independence. 

However, Spo modifies the evolution of the metric fluc- 
tuations according to the perturbed Einstein equation 
(|17p : and, in turn, this causes the corresponding evolu- 
tion of the matter perturbations through Eq. (|43l). We 
can write down the appropriate form of the perturbation 
equation as follows. First, we deflnc the "instantaneous" 
normalized densities at a cosmic time i, Qn = Pn/Pc 
where Pc = 3H^/8ttG is the critical density at the same 
instant of time. (Notice that the new densities equal the 
previously deflned ones, i.e. flj^ = i^N, only a,t t = to.) 
Making use of the deflnition of the relative DM fluctu- 
ation Sm and of the non-adiabatic perturbed pressure 
([39|) . we may cast Eq. (fT7|) in the following way: 



n H — n HmOm 

a a 



a 



(1 + ScDSd 



+ 9a^H{cl~cl) 



fc2 



(46) 



Next we use Eq. 
help of 

H'ia) = - 



to eliminate h from 



With the 



47rG 
aH{a) 
3H(a) 
2a 



[pM{a) + {1 + wja)) poia)] (47) 
We{a)r{a)' 



and 



51m (a) = 



l + r(a) ' 



1 + r(a) 



l + r{a) ' 



(48) 



(49) 



where r(a) is the "cosmic coincidence ratio" ((IJ between 
the DE and matter densities, we may flnally rewrite (I46p 
as follows 



1 



We (a) r{a) 
l + r{a) 



1 



5m (a) 



3 r(a) 



2 1 + r(a) 



(1 + 3ci) 



2sSDia) 



2 1 + r(a) a2 



+ 9Hia) (ci-ci) 



(a) 



fc2 



(50) 



If we would neglect the DE fluctuations {5d = 0, 0£) = 0), 
the r.h.s. of the previous equation would vanish. Un- 
der these conditions, one would be left with a decoupled. 
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second-order, homogeneous differential equation that de- 
termines the matter perturbations 5m- As could be ex- 
pected, the resulting equation coincides with Eq. (2.16) of 
Ref. [2^, where the approximation of neglecting the DE 
perturbations was made right from the start as an in- 
termediate procedure to investigate the amount of linear 
growth of the matter perturbations and put constraints 
on the parameter space of the AXCDM model. This pro- 
cedure was called "effective" in that reference, since all 
the information about the DE is exclusively encoded in 
the non-trivial EOS function — We{a). Let us write 
the homogeneous equation as follows, 



^ u 2 a 

(51) 

and let us assume a time interval not very large such that 
fi^) and We remain approximately constant. Looking for 
a power-law solution of (jSip in the limit 17 ^ 1, we 
find, for the growing mode. 



1-3(1-mJOd/5 



3(1 - We 



■ il n In ( 



(52) 

Since We < for any conceivable form of DE, this equa- 
tion shows very clearly that we should expect a growth 
suppression of matter perturbations (i.e. 6m ^ a", with 
n < 1) whenever a (positive) DE density Jlu is present 
within the horizon. Physically speaking, we associate this 
effect to the existence of negative pressure that produces 
cosmological repulsion of matter. 

However, being the DE density non-constant in general 
((5d 7^ 0), the DE perturbations themselves (and not only 
the value of the background DE density) should act as a 
source for the matter fluctuations. This effect is precisely 
encoded in the inhomogeneous part of Eq. (jSOp . i.e. in its 
r.h.s which is, in general, non-zero for 6d,Qd 7^ 0. In 
order to better appreciate this effect, let us consider an- 
other simplified situation where the analytical treatment 
is still possible, namely, let us assume an adiabatic regime 
(c^ = c^) with roughly constant EOS {we — const.) at 



very large scales (for which k in Eq. (j45p is very small, 
and hence the Oo component becomes negligible). Under 
these conditions, Eq. (|^^ greatly simplifies as follows: 



5'd = 



(1 



+ We) ^ _ Q 

aH 2 ^ ^ 



M J 



(53) 



where in the second step we have used . The rates of 
change of the perturbations for matter and DE, therefore, 
become proportional in this simplified setup. To be more 
precise, we see from ((53)) that for We > — 1 (quintessence- 
like behavior of the composite DE fluid) the matter fluc- 
tuations that are growing with the expansion > 0) 
trigger fluctuations in the DE also growing with the ex- 
pansion (S'jj > 0), whereas for < — 1 (phantom-like 
behavior of the DE) we meet exactly the opposite sit- 
uation, i.e., increasing fluctuations in the matter den- 
sity {S'jyj > 0) lead to decreasing fluctuations in the DE 



{S'jj < 0). Note that upon trivial integration of (|53|) one 
finds Sd — (1 + We) (5m + C', where C is a constant deter- 
mined by the initial conditions. For C — one obtains 
a result that fits with the well-known adiabatic initial 
condition relating the density contrasts of generic matter 
and DE components [1^ , 



WM 



WD 



(54) 



where, for non-relativistic matter, we have wm — 0, and 
lud is, in this case, the effective EOS We of the composite 
DE fluid. Since a positive DE density always leads to cos- 
mological repulsion, it follows from (|53p that one should 
expect some inhibition (resp. enhancement) of the mat- 
ter growth for the quintessence-like (resp. phantom-like) 
case. 

Although the previous example illustrates the impact 
of the DE fluctuations on the matter growth for a simple 
situation, a more complete treatment is required in the 
general case. In practice, this means that we have to nu- 
merically solve the system or, if desired, replace 
the last equation with the second order inhomogeneous 
Eq. ([50]) whose r.h.s. depends on the density contrast 
and the velocity gradient for the DE, 5d,0d 0. No- 
tice that the presence of overdensity DE perturbations 
{5d > 0) does not necessarily imply the inhibition of 
the corresponding matter perturbations since the coeffi- 
cient 1 -I- 3c^ in front of So on the r.h.s. of Eq. (|5D|) is 
positive for non-adiabatic DE perturbations. Only for 
= we meet the aforementioned possibility because 
~ We is usually negative, unless w^ is rapidly decreas- 
ing with the expansion, see Eq. ([^0]) . In this sense, the 
discussion above, based on Eq. ([55]) . is only valid at very 
large scales, specifically for fc-modes whose length scale 
I ^ is outside the sound horizon (cf. section |V|) ^. 
However, at smaller scales, specially at scales inside the 
sound horizon, and for a general non-adiabatic regime, 
we need to solve the aforesaid complete system of equa- 
tions for the basic set of perturbation variables for the 
metric, matter and DE: 



(/i(a, fc), SMia,k), Soia^k), 9D{.a,k)). 



(55) 



In this way, we have extended the effective treatment of 
the DE perturbations presented in Ref. [1^, and we are 
now ready to better assess the scope of its applicability. 
In section IVIIl we will apply this general formalism to 
the AXCDM model. 



^ In section Hxl we will see that this particular situation does ac- 
commodate very well the realistic dynamics of the cosmic per- 
turbations for matter and DE during the early stages of the evo- 
lution. 
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V. SOME GENERIC FEATURES OF THE DE 
PERTURBATIONS 

In the present section, we summarize some character- 
istic features of the DE perturbations. Many properties 
which are, in principle, common to any model with a 
self-conserved DE, will be later exemplified in section HXl 
within the non-trivial context of the so-called AXCDM 
model [13, il. 



A. Divergence at the CC boundary 

In general, the EOS of the DE will be a dynamical 
quantity, Wg = We{a)- In many models, the EOS may 
change from quintessence-like (—1 < We < ~l/3) to 
phantom (we < —1) behavior or vice versa, acquiring 
therefore the value = — 1 (also referred to as the 'CC 
boundary') at some instant of time. This is problematic 
since, as we shall see next, the equations for the pertur- 
bations diverge at that point. 

The divergence at the CC boundary is common to any 
DE model and has been thoroughly studied in the litera- 
ture (see e.g. [H, S^l). The problem can be readily seen 
by direct inspection of Eqs. ([331) and (I45|l . Note that, 
even though diverges at the crossing (cf. (1301)), the 
combination (1 + We) remains finite and, therefore, 
Ea. (|44ll is well-behaved. Thus, the problem lies exclu- 
sively in the (1 + We) factor in the denominator of (|^5l) . 
One might think that the divergence can be absorbed 
through a redefinition of the variables, but this is not the 
case. 

Getting around this difficulty is not always possible. It 
is well-known that there is no way for a single scalar field 
model to cross the CC boundarv [isj . The simplest way 
to avoid this problem is to assume two fields {Q,P), e.g. 
one that works as quintessence wq > — 1 and dominates 
the DE density until the CC-crossing point, and beyond 
it the other field retakes the evolution with a phantom 
behavior wp < —1, or the other way around; see [46, 47] 
for a detailed discussion and specific parameterizations of 
We- As we will see, in the AXCDM model the additional 
restrictions needed to avoid this divergence will further 
constrain the physical region of the parameter space. 



B. Unbounded growth for adiabatic DE 
perturbations 

Another well-known problem is the unbounded growth 
of the DE perturbations for a negative squared speed of 
sound cj.. As already mentioned in the previous section, 
in the adiabatic case we have ~ We, which is in general 
negative as long as the EOS parameter is not varying too 
fast. As a result, the adiabatic DE perturbations may 
lead to explosive growth unless extra degrees of freedom 
are assumed (see e.g. [l^l for discussion). 



In order to better see the origin of the problem, let us 
rewrite Eqs.(|331) ^^'^ SSI) in terms of conformal time rj, 
which is easily done by making use of / = H f (for 
any /): 




a (1 + We) 



As in section Hill we have defined Ti. — a/ a = Ha. If we 
use the two equations above and Eq. (jiSI to get a second 
order differential equation for 5b , we arrive at 



(58) 



where the second term on the r.h.s. represents other 
terms linear in these variables. Assuming that the var- 
ious perturbations are initially more or less of the same 
order, we see that the first term on the r.h.s. of ([551) '^i^l 
dominate provided 



Csdr] 



> 1 



(59) 



Notice that, for constant sound velocity, this condition 
simply tells us that the wavelength of the modes satisfies 
£ ^ ^ Xs, where As = CsTJ is the "sound horizon". 
Eq. (|59p is a generalization of this condition for arbitrary 
sound speed, in which case the sound horizon is given by 



A. 



Cgdrj 



Cg da 
a^H{d) 



(60) 



and constitutes a characteristic scale for the DE pertur- 
bations. As we will see next, the DE is expected to be 
smooth for scales well below it [i^ . 

For scales well inside the sound horizon, ([58l) becomes 
the equation of a simple harmonic oscillator, whose solu- 
tion is (in what follows, we assume constant for sim- 
plicity): 



5d = Cie 



iCs krj 



— icskri 



(61) 



where Ci and C2 are constants. We see that, for < 
(i.e. imaginary Cs) and neglecting the decaying mode, the 
DE perturbations grow exponentially. Obviously, this sit- 
uation is unacceptable for structure formation ^ . On the 
other hand, if > the DE density contrast oscillates. 
Since 6m grows typically as the scale factor a, this ensures 



® Let us note that this kind of problem need not occur when the DE 
is a pure running A , for then the perturbation Sp/^ is no longer 
an independent dynamical variable and is, in fact, determined 
by an algebraic function of the other perturbation variables, see 
\m for details. 
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that the ratio 5d/5m ''^ Sd/cl with the expansion. 
In other words, this tells us that the DE is going to be 
a smooth component (as it is usually assumed) as long 
as we are well inside the sound horizon. This feature is 
treated in more depth in the following section. 



Smoothness of DE below the sound horizon 



As a matter of fact, Eq. ([58|l is an oversimplification. 
In addition of having a term proportional to 6d, we also 
have one depending on its first derivative. So we can 
write that equation more precisely as follows: 



Sd = Di5d + D2 Sd + 0{h, 5m, Bd) 



(62) 



which gives us not just a simple, but a damped harmonic 
oscillator. The coefficients Di and D2 are, in general, 
functions of the conformal time 77. So the DE density 
contrast does not only oscillate, but its amplitude de- 
creases with time. Indeed, it was shown in [43,] that the 
quantity 



(63) 



which corresponds to the density contrast in the DE rest 
frame, oscillates with an amplitude A that decreases ac- 
cording to 



e)/2 



(64) 



The damped oscillations of the DE density contrast 
are clearly seen in the AXCDM model, as we will show 
in section Hxl 

Finally, we may ask ourselves whether the scales rel- 
evant to the LSS surveys [3] are well inside the sound 
horizon or not. Note that, in a matter dominated Uni- 
verse with negligible CC term and constant c^, we have 



simple form 



-"0 "m ' 



and the sound horizon (|60p takes the 



\s = 



2cs 



(65) 



Thus, in general, we expect that the size of the sound 
horizon at present (oq = 1) should be roughly of the or- 
der of the Hubble length H^^ . On the other hand, the 
observational data concerning the linear regime of the 
matter power spectrum lie in the range 0.01/iMpc"^ < 
k < 0.2hMpc~^ which corresponds to length scales 



e 



k comprised in the interval (600 i?o) 



< 



< 



{30 Hq)^^, hence well below the sound horizon (at least 
for not too close to 0). Therefore, according to the pre- 
vious discussions, we expect the DE density to be smooth 
at those scales, and indeed it will be so for the AXCDM 
model. Nevertheless, as we will see in section IIXI the 
larger the scale £ or the smaller the speed of sound Cg, 
the more important the DE perturbations are, because 
then (|551) is not such a good approximation. 



VI. THE AXCDM MODEL AS A CANDIDATE 
TO SOLVE THE COSMIC COINCIDENCE 
PROBLEM 

The AXCDM model [l3] provides an interesting way of 
explaining the so-called "cosmological coincidence pro- 
blem" (cf. section HI)) . The idea is related to the pos- 
sibility of having a dynamical component A, called the 
"cosmon" ^, which interacts with a running cosmological 
constant A. If the matter components are canonically 
conserved, the composite DE "fluid" made out of X and 
the running A will be a self-conserved medium too. The 
dynamics of the AXCDM universe is such that its com- 
posite DE may enforce the existence of a stopping point 
after many Hubble times of cosmological expansion. As a 
result, this modified FLRW-like universe can remain for a 
long while in a situation where the coincidence ratio ([T]) 
does not change substantially from the time when the 
DE became significant until the remote time in the fu- 
ture when the stopping point is attained. Subsequently, 
the Universe reverses its motion till the Big Crunch. 

The total DE density and pressure for the AXCDM 
universe are obtained from the sum of the respective CC 
and X components: 



Pd = PA + Px , 



Pd = PA+Px 



(66) 



The evolving CC density p\{t) = A(t)/87r G of the model 
is motivated by the quantum field theory formulation in 
curved space-time by which the CC is a solution of a 
renormalization group equation. Following [l^, [l^, [l^, 
[H,!!^, the CC density emerges in general as a quadratic 
function of the expansion rate: 



Pa{H) 



g Ml {H' Hi) , (67) 



where = pa{H = Hq) is the present value. The di- 
mensionless parameter ly is given by 



12 n Ml 



(68) 



where M is an effective mass parameter representing the 
average mass of the heavy particles of the Grand Unified 
Theory (GUT) near the Planck scale, after taking into 
account their multiplicities. Depending on whether they 
are bosons, or fermions, cr = -t-l,orcr = — 1, respectively. 
For example, for M = Mp one has \v\ =1^0, where 



uo = ^ '^2.Qx 10- 

iZTT 



(69) 



Originally, the cosmon appeared in Ref. [31 as a scalar field linked 
to the mechanism of dynamical adjustment of the CC. In the 
present context, the entity X is also differentiated from the CC, 
but if taken together they form a composite and interactive DE 
medium. 
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On physical grounds, we expect that this value of 
should be the upper bound for this parameter. In the 
next section, we will see if we can pinpoint a region of 
parameter space compatible with this expectation. 

The energy density associated to the cosmon compo- 
nent X is obtained from the total DE conservation law 
([36)1 and the composite form ((66)) . 



p'xia)+p'A{a)^ — {l + wx)pxia), (70) 
a 

where wx is the effective EOS parameter of X, 



Px = wxpx 



(71) 



In principle, wx could be a function of the scale factor. 
However, a simpler assumption that allows us to perform 
a completely analytic treatment, is to consider that the X 
component behaves as a barotropic fluid with a constant 
EOS parameter in one of the following two ranges: lox ^ 
— 1 (quintessence-like cosmon) oi ujx ^ — 1 (phantom-like 
cosmon). On the other hand, the EOS parameter for the 
running A component still remains as the cosmological 
constant one, wa = — 1, i.e. 



PA 



'PA ■ 



(72) 



From these assumptions, it is easy to find the following 
relation between the effective EOS parameter of the total 
DE fluid (|57)) and the EOS parameter of the cosmon wx- 



(1 + We{a)) pnia) = (1 + wx) Px{a) ■ 



(73) 



The normalized density of the cosmon component, 
flxio,) — Px(i)/p2, can be obtained from the previous 
relations after solving the differential equation 170]) . In 
this equation, we have p'j^{a) = (3i^/87r) Mp dH'^/da from 
(IS7)) . and the derivative dH^/da — 2H{a)H'{a) can be 
exphcitly computed from ([iT)) upon using (|75|) and 
One flnally obtains the differential equation 

^'xia) + - (l+ujx-e) Qxia) ^Si^nlja-^ . (74) 
a 

With the boundary condition that the current value of px 
is , the solution of ([7^ can be written in the following 
way: ® 



nx{a) 



(1 + 6) a-3(i+'"-^-^) -fofl- 



(75) 



where in the above equations we have used the notations 

(76) 



v{l + Wx) ■ 



(77) 



As we will discuss in more detail below, the parameter 
e must remain small (|e| < 0.1) in order to be compati- 
ble with primordial nucleosynthesis. For 1^ = the CC 
density ((67|) becomes constant. In this case, the two pa- 
rameters (|76|) and ((77|) vanish and Eq. (|75p boils down 
to the simplest possible form, which is characteristic of a 
self-conserved monocomponent system. 



Qxia) = ^ 



„-3{l+wx) 
Xtt 



(78) 



It is only in this particular situation where the cosmon 
X could be a self-conserved scalar field with its own dy- 
namics. But in general this is not so because in QFT 
in curved space-time we have good reasons to expect a 
running pA [H, [H, [III, and hence v ^ 0. Therefore, if 
the total DE is to be conserved, the dynamics of X is not 
free anymore and becomes determined as in ([75]) . 

The normalized total DE density fli) = po / p^c given 

by 



wxnl 



M 



V Wx 
— V 



-3{l+wx-f-) 



1 



Wx 



(79) 



where the various current normalized densities satisfy the 
relation VPf^j + = + rPf^ + n\ = I, which may be 
called the "AXCDM cosmic sum rule". Using ([75)) . the 
effective EOS of the DE in the AXCDM model can now 
be obtained explicitly. 



u;e(a) = -l + (l+c.^)§44' 



(80) 



with Vlx{a) and ^Id{o) given by ([75|) and (fTQ]) . respecti- 
vely 

The total DE density ((79)) varies in such way that the 
ratio ll]) can remain under control, which is the clue for 
solving the coincidence problem in a dynamical way [l3] ■ 
Indeed, the explicit computation of such ratio yields 



r{a) 



^li (1 - ^) 



Wx 

Wx - e 
e 

I 

Wx - 



-3(t«x-e) 



(81) 



The fact that the evolution of the cosmon X is completely deter- 
mined by the dynamics of the running pA (I67I I. together with the 
hypothesis of total DE conservation I I70I I. implies that X cannot 
be generally assimilated to a scalar field, which has its own dy- 
namics. In fact, as we have already mentioned, X is to be viewed 
in general as an effective entity within the context of the effective 
action of QFT in curved space-time. 



and it can be bounded due to the existence of a maximum 
(triggered by the ~ a~3(iux-c) term in the previous for- 
mula). Moreover, r{a) given above stays relatively con- 
stant (typically not varying more than one order of mag- 
nitude) for a large fraction of the history of the Universe 
and for a significant region of the parameter space p^j . 
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In contrast, in the standard concordance ACDM model, 
the CC density remains constant, pA = Pa' ^^'^ 
incidence ratio grows unstoppably with the cubic power 
of the scale factor, r{a) — H.'j^ / . In this scenario, 
it is difficult to explain why the constant pA = Pa 
the same order of magnitude as the matter density right 
now: p^j. Let us point out that the standard model ratio 
is just that particular case of ([5T|) for which i/ = (no 
running CC) and fl^ = (no cosmon). 

Before closing this section, we would like to make a re- 
mark and some discussion concerning the quadratic evo- 
lution law ([S7)) for the cosmological term. This equa- 
tion was originally motivated within the framework of 
the renormalization group (RG) of QFT in curved space- 
time [H, [H, [3, m, mi (see also fld\ for a short review). 
We point out a criticism against this approach that re- 
cently appeared in the literature [s^ . While the nature 
of this criticism was amply rejected in [l^ (see below for 
a summary) , it is fair to say that the question of whether 
a rigorous RG approach in cosmology is feasible is still 
an open question and remains a part of the CC problem 
itself. Although it is not our main aim here to focus 
on this fundamental issue, let us briefly sketch the situ- 
ation along the lines of Ref. [l^ , to which we also refer 
the reader for a summary of the rich literature proposing 
different RG formulations of cosmology both in QFT in 
curved space-time and in Quantum Gravity. 

The RG method in cosmology treats the vacuum en- 
ergy density as a running parameter and aims at finding a 
fundamental differential relation (renormalization group 
equation) of the form 



dpA 
din ji 



(82) 



which is supposed to describe the leading quantum con- 
tributions to it, where /3a is a function of the parameters 
P of the effective action (EA) and /i is a dimensional 
scale. The appearance of this arbitrary mass scale is char- 
acteristic of the renormalization procedure in QFT owing 
to the intrinsic breaking of scale invariance by quantum 
effects. The quantity pA in (|82p is a (/i-dependent) renor- 
malized part of the complete QFT structure of the vac- 
uum energy. Depending on the renormalization scheme, 
the scale /x can have a more or less transparent physi- 
cal meaning, but the physics should be completely in- 
dependent of it. Such (overall) /i- independence of the 
observable quantities is actually the main message of the 
RG; but, remarkably enough, the /^-dependence of the 
individual parts is also the clue of the RG technique to 
uncover the leading quantum effects. 

Essential for the RG method in cosmology is to un- 
derstand that, in order for the vacuum energy to acquire 
dynamical properties, we need a nontrivial external met- 
ric background. The dynamical properties of this curved 
background (e.g. the expanding FLRW space-time, char- 
acterized by the expansion rate H) are expected to induce 
a functional dependence pA — Pa{H)- The latter should 
follow from parameterizing the quantum effects with the 



help of the scale pL and then using some appropriate cor- 
respondence of p with a physical quantity, typically with 
H in the cosmological context, although there are other 
possibilities [ll,[ll [13. In this way, one expects to esti- 
mate the subset of quantum effects reflecting the dynam- 
ical properties of the non-trivial background. Although 
p cancels in the full EA, the RG method enables one to 
separate the relevant class of quantum effects responsible 
for the running. The procedure is similar to the RG in 
a scattering process in QCD; the parametrization of the 
quantum effects in terms of p is the crucial strategy to 
finally link them with the energy of the process through 
the correspondence p q (where g is a typical momen- 
tum of the scattering process) at high energy. One can 
also proceed in the same way in QED and electroweak 
theory (although here one can adopt more physical sub- 
traction schemes, if desired). The RG technique can ac- 
tually be extended to the whole Particle Physics domain. 
In cosmology, however, the situation is more complicated, 
partly because (as remarked above) the physical scale be- 
hind the quantum effects is not obvious. Still, one expects 
that it should be related with the expanding metric back- 
ground and hence the expansion rate H can be regarded 
as a reasonable possibility. On this basis, the heuristic 
arguments exhibited in [ll [ll, [H [H, [Uliil combined 
with the general covariance of the EA suggest that the 
solution of the RG equation (|82p should lead to the kind 
of quadratic law (|67p that we have used. 

According to [lB| , the point of view of Ref. [lOl is in- 
correct on two main accounts: first, because they try 
to disprove the running through the overall cancelation 
of the arbitrary scale p in the EA; and second, because 
they neglect the essential role played by the non-trivial 
metric background. As emphasized in [l^, the cancela- 
tion of p in the EA cannot be argued as a valid criticism 
because this fact is a built-in feature of the RG and it 
was never questioned. If this would be a real criticism, it 
would also apply to QED, QCD or any other renormaliz- 
able QFT, and nevertheless this is no obstacle for using 
the RG method in these theories as an extremely useful 
strategy to extract the dependence of the quantum effects 
on the physical energy scale of the processes, in partic- 
ular the so-called running coupling constants gs = gsil) 
and e = e{q) of the strong and electromagnetic interac- 
tions. Moreover, in the absence of a non-trivial metric 
background, there is no physical running of the vacuum 
energy, even though there is still /i-dependence of the 
vario us p arts of the EA and in particular of the CC, see 
e.g. [5l|. Therefore, at the end of the day such criticism 
seems to go against the essence of the RG method and 
its recognized ability to encapsulate the leading quantum 
effects on the physical observables. 

In cosmology, the principles of the RG should be the 
same, but there are two main stumbling blocks that 
prevent from straightforwardly extending the method in 
practice [iBl , to wit: i) the aforesaid lack of an obvi- 
ous/unique correspondence of p with a cosmological scale 
defining the physical running, and also (no less impor- 
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tant) ii) the huge technical problems related with the 
application of the RG within a physical (momentum- 
dependent) renormalization scheme in a curved back- 
ground. These difficulties are unavoidable here because 
we are dealing with QFT in the infrared regime and more- 
over the metric expansions cannot be performed on a flat 
background; indeed, there cannot be a flat background 
in the presence of a cosmological term! While these two 
problems remain unresolved in a completely satisfactory 
manner, it is legitimate to use the phenomenological ap- 
proach and the educated guess (e.g. the general covari- 
ance of the EA) to hint at the running law. This is the 
guiding principle followed in the aforementioned refer- 
ences and that led to Eq. ((GT]) . 

Finally, let us emphasize that irrespective of whether 
such law can be substantiated within the strict frame- 
work of the RG, the present study remains perfectly use- 
ful simply treating (j67p as an acceptable type of a phe- 
nomenological variation law and keeping also in mind 
that adding the cosmon may contribute to the resolution 
of the pressing cosmic coincidence problem. 



VII. DARK ENERGY PERTURBATIONS IN 
THE AXCDM MODEL 

In this section, we further elaborate on the conditions 
to bound the ratio r = r(a) and discuss the constraints 
on the parameter space, in particular the impact of the 
DE perturbations on these constraints. In section |Vl we 
discussed analytically some generic features about DE 
perturbations. In principle, those results should apply 
to any model in which the DE is self-conserved. The 
AXCDM model, given its peculiarities (a composite DE 
which results in a complicated evolution of the effective 
EOS), constitutes a non-trivial example of that kind of 
models. In this sense, it is interesting to use the AXCDM 
to put our general predictions to the test. At the same 
time, this will allow us to impose new constraints on the 
parameter space of the model, improving its predictivity. 

The parameter space of the AXCDM model was al- 
ready tightly constrained in [2^ . In that work, the mat- 
ter density fluctuations were analyzed under the assump- 
tion that the DE perturbations could be neglected. As 
a first approximation this is reasonable since, as we have 
discussed in section |Vl the DE is expected to be smooth 
at scales well below the sound horizon. Thus, we will 
take the results of [1^ as our starting point and will 
check numerically the goodness of that approximation. 
Finally, we will further constrain the parameter space 
using the full approach presented in the present work. 
Let us summarize the constraints that were imposed in 

1. Nucleosynthesis bounds: As already commented, 
the ratio (|8ip between DE and matter densities 
should remain relatively small at the nucleosynthe- 
sis time, in order not to spoil the Big Bang model 



predictions on light-element abundances. Requir- 
ing that ratio to be less than 10%, it roughly trans- 
lates into the condition |e| < 0.1, where the param- 
eter e was defined in Eq. ([77|) . 

2. Solution of the coincidence problem: In Ref. [l3|, 
where the AXCDM model was originally introduced 
as a possible solution to the coincidence problem, 
it was shown that there is a large sub-volume of 
the total AXCDM parameter space for which the 
ratio r(a) remains bounded and near the current 
value To (say, \r{a)\ < 10 tq, where tq ~ 7/3) dur- 
ing a large fraction of the history of the Universe. 
Thus, the fact that the matter and DE densities 
are comparable right now may no longer be seen 
as a coincidence. Such solution of the coincidence 
problem is related to the existence of a future stop- 
ping (and subsequent reversal) of the Universe ex- 
pansion within the relevant region of the parameter 
space. 

3. Current value of the EOS parameter: Recent stu- 
dies (see eg. (H) suggest that the value of the DE 
effective EOS should not be very far from —1 at 
present. Although these results usually rely on the 
assumption of a constant EOS parameter (and thus 
are not directly applicable to the AXCDM model), 
we adopted a conservative point of view and stuck 
to them by enforcing the condition \1 + We{a = 
1)1 < 0.3 on the EOS function ([50]) . 

4. Consistency with LSS data: As said before, in 
(23} we studied the growth of matter density fluc- 
tuations under the assumption that the DE was 
smooth on the scales relevant to the linear part of 
the matter power spectrum. From the fact that the 
standard ACDM model provides a good fit to the 
observational data, we took it as a reference and 
imposed that the amount of growth (specifically, 
the matter power spectrum), of our model did not 
deviate by more than 10% from the ACDM value 
("F-test" condition). This condition can be also 
be justified from the observed galaxy fluctuation 
power spectrum, see [IJ for more details. 

The upshot of that analysis was that there is still a big 
sub- volume of the three dimensional AXCDM parameter 
space {v, uJx, satisfying simultaneously the above 
conditions^. The projections of that volume onto the 
three perpendicular planes {ly, {v, ujx) and (fl^,LUx) 
are displayed in Figs.[T]and[2] (shaded regions). These re- 
gions where already determined in Ref. 23]. In the next 
section, we will discuss how the final set of allowed points 



^ In Ref. [2^ , we took a prior for the normalized matter density at 
present, specifically Q,^ = 0.3. This means that H'^ = 0.7 for a 
spatially fiat Universe. For better comparison with those results, 
we keep this prior also in the present work. 
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becomes further reduced when we take into account the 
analysis of the DE perturbations. 



A. Divergent behavior at the CC boundary 

As discussed in section |Vl if the effective EOS of the 
model crosses the CC boundary {we = —1) at some point 
in the past, the perturbation equations will present a 
real divergence. Obviously, this circumstance makes the 
numerical analysis unfeasible at the points of parameter 
space affected by the singularity. 

In the absence of an apparent mechanism to get around 
this singularity, we are forced to restrict our parameter 
space to the subregion where the solution of the perturba- 
tion equations (I43p - (l46p is regular, namely by removing 
those points of the parameter space that present such a 
crossing in the past, because these points can not be- 
long to a well defined history of the Universe. In the 
absence of a more detailed definition of the cosmon en- 
tity X, this new constraint is unavoidable. This should 
not be considered as a drawback of the model, for even in 
the case when one uses a collection of elementary scalar 
fields to represent the DE, one generally meets the same 
kind of divergent behavior as soon as the CC boundary 
is crossed, unless some special conditions are arranged. 
In other words, even if the components of the DE are as 
simple as, say, elementary scalar fields with smooth beha- 
vior and well-defined dynamical properties (including an 
appropriately chosen potential), there is no a priori guar- 
antee that the CC boundary can be crossed safely 45|. 
It is possible to concoct ingenious recipes, see e.g. 46| . 
such that the perturbation equations become regular at 
the CC boundary, but the procedure is artificial in that 
one must introduce new fields (one quintessence-like and 
another phantom-like) satisfying special properties such 
that their respective EOS behaviors match up continu- 
ously at the CC-crossing. Apart from the fact that fields 
with negative kinetic terms are not very welcome in QFT, 
one cannot just replace the original fields with the new 
ones without at the same time changing the original DE 
model! As we will see below, in the AXCDM case the 
absence of CC-crossing projects out a region of the pa- 
rameter space which is significantly more reduced, and 
therefore the predictive power of the model becomes sub- 
stantially enhanced. 

In section 6 of the first reference in [l3| , it was shown 
that the necessary and sufficient condition for having a 
CC boundary crossing in the past within the AXCDM 
model is that the parameter b given in (j76p is positive. As 
can be readily seen, this will happen whenever v and fl'j^ 
have the same sign (where we use the fact that ujx < 
and |e| ^ \oJx\ in the relevant region of parameter space). 
From the cosmic sum rule of the AXCDM model, we have 
= 1 — r2"j — r2° ; thus, using our prior fll^ — 0.3, 
the set of allowed points (for which fl^ has different sign 
from v) are those comprised in the striped areas in Figs.[T] 
and [21 This leaves us with a very small allowed region in 




FIG. 1: Projection of the 3D physical volume of the AXCDM 
model onto the — 0° plane. Shaded area: points that satisfy 
all the constraints in |231 , see also section IVIII of the present 
work. Striped area: points that are not affected by the diver- 
gence at the CC boundary discussed in Sect. I VII Al The final 
allowed region is the one both shaded and stripped. As a re- 
sult of considering the DE perturbations, the possible values 
of the parameters become strongly restricted, which implies a 
substantial improvement in the predictive power of the model. 



each plane, which is just the corresponding intersection 
of the shaded area and the striped one. 

At the end of the day, it turns out that most of the 
points in the shaded area in Figs.[l] and [2] (viz. those 
allowed by the conditions stated in the previous section 
and the analysis of [1^) are ruled out by the new con- 
straint emerging from the DE perturbations analysis, and 
hence we end up with a rather definite prediction for the 
values of the AXCDM parameters. In particular, we find 
from these figures that only small positive values of v are 
allowed, at most of order v ^ 10^^. Let us emphasize 
that this is in very good agreement with the theoretical 
expectations mentioned in section IVTl Recall that, from 
the point of view of the physical interpretation of v in 
Eq. we expected v in the ballpark of i^o 10~^ at 
most - see Eq. ([69]) - since the masses of the particles 
contributing to the running of the CC should naturally 
lie below the Planck scale Let us mention that the 
interesting bounds on v obtained in Ref. [H^ on the basis 
of the so-called generalized Second Law of gravitational 



Let us clarify that the tighter bounds on u determined in Ref. [4l| 
are possible only because, in the latter work, the DE is not con- 
served and there is no cosmon. As we have shown in l23t , a run- 
ning cosmological constant model without a self-conserved DE 
cannot solve the coincidence problem in a natural way because 
the required values of v are too large and, hence, incompatible 
with the physical interpretation of this parameter. 
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FIG. 2: Projection of the 3D physical volume of the AXCDM 
model onto the u — ujx and 0° ~ planes. Shaded: points 
that satisfy all the constraints in 23]. Striped: points that, 
in addition, are not affected by the divergence at the CC 
boundary discussed in Sect. IVII Al The final allowed region 
is the one both shaded and stripped. 



thermodynamics would suggest that only the eflFective 
mass near the Planck mass is allowed. However, let us 
point out that such study has been performed without 
including the non-trivial effect from the cosmon. 

A very important consequence of the dark energy per- 
turbative constraint is that the effective EOS of the DE 
can be quintessence-like only, i.e. —1 < We < —1/3. To 
prove this statement, let us start from Eq. ([75]) . For the 
current values of the parameters, this equation can be 
rewritten as 



first, the cosmon component is necessarily phantom-like 
{u!x < —1) in the allowed region by the DE perturba- 
tions; and, second, its energy density at present is nega- 
tive; namely, = 0.7 — 51° < because from Fig.[2h> 
we have 57^ > 0.7. Therefore, since the r.h.s. of ((83)) is 
constrained to be positive and il^j = 0.7 > 0, we are en- 
forced to have wj? > —1. However, the fulfilment of this 
condition at present implies its accomplishment in the 
past, i.e. We{a) > —1 (Va ^ 1), otherwise there would 
have been a crossing of the CC boundary at some earlier 
time, which is excluded by the analysis of the DE per- 
turbations. The upshot is that the EOS of the DE in 
the AXCDM model can only appear effectively as quin- 
tessence (q.e.d). In reality, it only mimics quintessence, 
of course, as its ultimate nature is not; such DE medium 
is a mixture made out of running vacuum energy and a 
compensating entity that insures full energy conservation 
of the compound system. 

Dark energy components X with negative energy den- 
sity are peculiar in cosmology since, in contrast to stan- 
dard DE components, they satisfy the strong energy con- 
dition (like ordinary matter), and as a result the gravita- 
tional behavior of X is attractive rather than repulsive. 
Due to this double resemblance with matter and phan- 
tom DE (although with the distinctive feature px < 0), 
such components can be called "phantom matter" [17| . 
Being X in general an effective entity, such "phantom 
matter" behavior is actually non-fundamental. 



B. Adiabatic speed of sound in the AXCDM 

In the equations we assumed the most gen- 

eral case in which the perturbations could be non- 
adiabatic. Moreover, we have shown that the adiabatic 
case usually leads to an unphysical exponential growth 
of the perturbations as a result of in being neg- 
ative. Next we will check that, indeed, the most com- 
mon situation in the AXCDM model is to have < 0. 
Notwithstanding, adiabatic perturbations are not com- 
pletely forbidden in the present framework, as there is a 
small region of the parameter space for which could 
be positive. 

From Eqs. (|^ and ([75)) . and making use of the DE 
conservation law p6p . we find that the adiabatic speed 
of sound for the AXCDM model can be cast as follows, 



3 Six (a) 



(84) 



With the help of (|75|) . we can rewrite the last expression 
as 



1 



n% = (1 



-u;x^ 



(83) 



= w^la = 1) is the value of the effective 
EOS parameter at the present time. Looking at Figs.[2^ 
and ^p, we realize the following two relevant features: 



clia) 



fix (a 



■(1 + h){ujx - e) a-3(i+-A-.) (85) 



We want to find out the condition for this expression to 
be positive. As {lux — e) < (remember that ujx < 
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and |e| < 0.1, due to nucleosynthesis constraints), that 
condition simply reads 



fix (a) 



(1 + fe) < 0. 



(86) 



The cosmon energy density fix (a) cannot vanish be- 
cause, in such case, the perturbation equations would di- 
verge. Indeed, lijc (a) = corresponds to a CC-boundary 
crossing at some value of the scale factor in the past, cf. 
Eg. ([M)) . Thus, being ilx{a) a continuous function, it 
must have the same sign in the past as that of its present 
value, i.e. il^ /flx{a) > (Va ^ 1)). In short, the final 
condition that ensures that > is 



(l + b) <0. 



(87) 



If this condition would be satisfied, then > would 
hold for the entire past history of the Universe and, im- 
der these circumstances, the adiabatic equations may be 
used. It turns out that the relation ([57|) can be satisfied 
in the AXCDM model, although only in a narrow range 
of the parameter space. In fact, from the definition of the 
parameter b in ([7^ . the expectation that \llJx\ — 0(1), 
and neglecting e, we see that (j87l) is approximately equiv- 
alent to 



^0^ > 



(88) 



Given the fact that v was found to be positive and small 
(cf. Fig.[l][2]) and, at the same time, Jl^f < (see previ- 
ous section), the above condition does not leave much 
freedom within the allowed parameter space, roughly 
—V < Vl\ < 0. This narrow strip is, however, not neces- 
sarily negligible; e.g. if we take v of order of vq ~ 10~^ 
(cf. Eq. ([69]) ). this possibility is still permitted in the 
parameter space, see Figs.[l][2] In such case, the present 
cosmon density could still be of the order or larger (in ab- 
solute value) than, say, the current neutrino contribution 
to the energy density of the Universe (0° ~ 10""^). No 
matter how tiny is (in absolute value) a negative cosmon 
contribution to the energy density, it suffices to take care 
of the cosmic coincidence problem along the lines that 
we have explained. Therefore, the adiabatic contribu- 
tion is perfectly tenable, but the numerical analysis of 
the subsequent sections remains essentially the same (as 
we have checked) independently of whether the sound 
speed of the DE medium is adiabatic or not. For this 
reason, in what follows we will assume the more general 
situation of non- adiabatic perturbations, with the under- 
standing that adiabatic ones can do a similar job in the 
corresponding region of the parameter space. 



VIII. THE MATTER POWER SPECTRUM 

In this section, we compare the matter power spec- 
trum predicted by the AXCDM model with the observed 



galaxy power spectrum measured by the 2dFGRS sur- 
vey The AXCDM matter power spectrum is found 
by evolving the perturbation equations from 
a — a.i to the present (oq = 1), where <C 1 is the scale 
factor at some early time, but well after recombination. 
In these equations, we must of course use the expansion 
rate ^ with the fuU DE density dH]). 

In order to set the initial conditions at a = a^, we use 
the prediction from the standard ACDM model. Indeed, 
the standard ACDM model provides a good analytical fit 
to the 2dFGRS observed galaxy power spectrum. Taking 
this fit as our starting point, we compute analytically the 
values of the ACDM perturbations at an arbitrary scale 
factor. Since the DE does not play an important role 
until very recently, we may assume that the initial mat- 
ter and metric perturbations at a = for the AXCDM 
model are the same as for the ACDM model. 



A. Initial matter and metric perturbations 

As previously commented, the perturbed equations are 
evolved from some ^ 1 to aq = 1 . The value of ai is 
unimportant, provided that it lies well after recombina- 
tion (to insure that all the processes encoded in the trans- 
fer function have already taken place). For definiteness, 
we take — 1/500, i.e. cosmological redshift z ~ 500. 
Next we specify the initial conditions at a = o^. 

For the matter and metric perturbations, the initial 
conditions in the standard cosmological model can be 
computed analytically [s^ [s^ . The matter perturbed 
equations in the standard ACDM model are and 
([50)1 taking ^^(a) = Od{o) = 0. In these conditions, the 
r.h.s. of the differential equation ((50|) vanishes and we 
are left with an homogeneous equation. In the standard 
model, this equation is just Eq. (I5ip with We = —1, i.e. 



6l,{a) + l{2-nM{a))^-^ 



3 ^ , .hiia) 



= 0. 
(89) 

Since matter is conserved, we have r2jv/(a) — 
{Ho/H{a)fn%/a^. Moreover, from Eq. (again 
with We = —1) one can also see that Q.Mio) — 
— {2a/3) H'{a)/H{a). Hence, the homogeneous equation 
([89]) can be conveniently cast as follows: 



(90) 

This equation is now in a standard form and can be 
solved analytically [1^ . Let us first introduce the variable 
D{a) = 6M{a)/Sroi (the so-called growth factor), where 
i5ref is the matter density contrast at some initial scale 
factor. In the initial matter era {a — a.i <ti I) wherein 
the cosmological term is negligible, we have H'^{a)/H'^ — 
very good approximation; then, by imposing 
the boundary conditions D(a) cx a and D'{a) — 1, the 
growing solution of (|90p is simply D{a) = a, as can be 
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easily checked. In the general case, one can find the solu- 
tion for the growing mode which reduces to the previous 
one deep into the matter dominated era. The final result 
reads 



D{a) = 



(5m (a) _ 5r!°,i7(a) T 



3rcf 



da 



Ho Jo {hH{h)/Hof 



(91) 



In particular, for H{a)/Ho — (il"^)-^/^ a~'^/^ it just boils 
down to the solution D{a) = a corresponding to the early 
matter dominated epoch, as expected. 
Furthermore, Eq. (l43l) leads to 



(92) 



h{a) = 2aH{a)^^5M{a), 
D{a) 

and for D{a) = a it renders the initial condition 
h{ai) = 2H{ai)5M{ai) 



(93) 



for the metric fluctuation. Later on, when the DE (i.e. 

> in the ACDM) starts to play a role, the mat- 
ter (and metric) fluctuations become suppressed. The 
suppression is given by the value of the growth factor 
D{a), which is no longer proportional to the scale fac- 
tor From ([?T|) it is clear that 5m{o)/ D{a) is a con- 
stant, which can be written as SM{o-i)/ai at early times 
(when DlcLi) = Ui) and as 5A/(ao)/I?(ao) at the present 
time. Therefore, 

SAiiai) ^ -^SAiiao) , (94) 

h{ai) ^2H^{ai)^5M{ao), (95) 

where Do = D{ao) and the subindex in H\(ai) has 
been added to emphasize that the initial value of the 
Hubble parameter is to be computed within the ACDM 
model. Note that, instead of setting the value of 
h(ai), we could have chosen to put initial conditions 
on the derivative of the density contrast, S'j^j. In that 
case, as it is evident from (|M)) . we would have that 
^'ivii'^i) ~ ^m{olo) I Do- Then the initial value of the 
metric fluctuation is constrained by Eq. to be 

h{ai) = 2Hj^x{ai)aiSM{ao)/Do, where now HAxiat) is 
the AXCDM value of the Hubble function. Note that this 
value of h(ai) is not exactly the same as that in (j95p . since 
H\{ai) and i?Ax(ai) are not identical. However, being 
the difference rather small (as we have checked numeri- 
cally), the behavior of the perturbations does not depend 
signiflcantly on that choice. 



Notice that if Ha is small and essentially constant, then the 
growth factor D{a) takes the approximate form D{a) a", with 
n = 1 — GOa/S < Ij S'S it follows from II52I I for = —1, or 
from H89|l . This demonstrates, if only roughly, the suppression 
behavior in an explicit analytic way. In the general ACDM case, 
however, the solution for the growing mode is given by ((91} , in 
which H is the full expansion rate of the standard model. 



The equations ((94|) and ([95]) give us the initial condi- 
tions at ai — 1/500 for the matter and metric perturba- 
tions in terms of the density contrast today (5m (ao). We 
associate the latter with the 2dFGRS observed galaxy 
power spectrum fitted in the ACDM model, as detailed 
below. 

The matter power spectrum of the ACDM model can 
be approximated as |41i] : 



PA(fc) = |(5M(fc)P ^ Ak T\k) 



(96) 



where ftj^ 



"a/ 



ri^. It assumes a scale-invariant 



(Harrison-Zeldovich) primordial spectrum, as generically 
predicted by inflation. This primordial spectrum is mod- 
ifled when taking into account the physical properties of 
different constituents of the Universe, in particular the 
interactions between them. All these effects are encoded 
into the scale-dependent transfer function T{k), which 
describes the evolution of the perturbations through the 
epochs of horizon crossing and radiation/matter transi- 
tion. The growth at late times which, in the ACDM 
model, is independent of the wavenumber, is described 
by the growth function g{^). Finally, y4 is a normaliza- 
tion factor. 

The transfer function can be accurately computed by 
solving the coupled system formed by the Einstein and 
the Boltzmann equations. Although a variety of numeri- 
cal fits have been proposed in the literature, here we use 
the so-called BBKS transfer function ^54| |: 



T(fc) = 



ln(l + 2.34g) r 



where 



2Mq 
(5.46g)3 

l = q{k) = 



-1/4 



I + 3.89q + {m.lqY 
(6.71g)4T 

k 



(hT) Mpc"^ 
and F is the Sugiyama's shape parameter [SE 



(97) 



(98) 



(99) 



On the other hand, for the growth function we assume 
the following approximation [56] : 



^)(i + ^ 

2/V 70 



(100) 

which reflects the suppression in the growth of per- 
turbations caused by a positive cosmological constant. 
The normalization coefficient A is related to the CMB 
anisotropics through [4l| 



67r2 Q 



A = i2iHr ^ 



rjns — PS 

-^0 



(101) 
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where Qrms-ps is the quadmpole amplitude of the CMB 
anisotropy (see below for more detailed explanations), 
Ih = H-^ ~ 3000/i^iMpc is the Hubble radius and Tq ~ 
2.725K, the present CMB temperature. Therefore, the 
value of the normalization factor A could in principle 
be inferred from measurements of the CMB. However, 
we have obtained it by fitting the power spectrum (|96p 
to the 2dFGRS observed galaxy power spectrum 4], as 
discussed below. 

We assume h = 0.7 for the reduced Hubble parameter 
and a spatially flat Universe with fi^j — 0.3 (hence — 
0.7 for the flat ACDM model) in order to be consistent 
with our assumption in previous analyses p7ll23|. The fit 
to the 2dFGRS observed galaxy power spectrum P2dF{k) 
is obtained assuming the matter budget composed of a 
baryonic part fl^ = 0.04 and a dark matter contribution 
^DM = 0.26. In order to calculate the best fit, we use 
the formula (|96p to minimize the ^-square distribution 

P.WP 

in terms of the normalization A. There are 39 values of k 
in the 2dFGRS data, so the number of degrees of freedom 
is ndof — 38. We find, as best fit, the value 

A = 8.99 X 10^/i-'*Mpc^ (103) 

with = 0.43. From (|10ip we see that this value 
of A implies Qrms~ps — 20.85 fiK. Let us now 
clarify that Qrms-PS is not the observed quadrupole 
CMB anisotropy (usually denoted Qrms), but rather 
the value derived from a fit to the entire CMB power 
spectrum (PS). For a power-law spectrum with n = 1 
(i.e. for a scale- invariant PS), the COBE team ob- 
tained Qrms-PS = 18 ± 1.6 fiK [st'I, and moreover they 
found that the observed Qrms is smaller than the fitted 
Qrms-PS- Whether this is a chance result of cosmic vari- 
ance or reflects the physical cosmology is not known [s^ . 
Our fitted value for Qrms-PS isXls within the 2a range 
of the corresponding COBE value (although when the 
quadrupole itself is not used in the fit, the COBE un- 
certainties become larger (l^l). As several authors have 
noted [55I [ssl , such a normalization may be inadequate 
for models with a cosmological constant, given the fact 
that the CMB spectra of A-dominated models is quite dif- 
ferent from a simple power-law, specially at large scales 
(low multipoles). For instance, in \58\ it is proposed an 
alternative normalization for the ACDM model, which 
for h = 0.8 and = 0.7 yields Qrms-PS = 22.04 ^iK, 
with an error of the order of 11%, which is in agreement 
with our result. In general one can find a number of dif- 
ferent values for Qrms-PS in the literature depending on 
the kind of analysis performed or the data set used, and 
this is why we preferred to compute the normalization di- 
rectly from a fit to the matter power spectrum. Finally, 
let us emphasize that our value for Qrms-PS lies within 
the 95% confidence interval for the observed quadrupole 
anisotropy {Qrms) by both COBE |57|] and WMAP 



Therefore, we will assume (l94l) and (|95)) as the initial 
conditions for the matter and metric perturbations, iden- 
tifying (SA/(ao) with Suik) from the formula using 
the fitted coefficient p03p . 

B. The AXCDM matter power spectrum 

The procedure discussed above helps us to set initial 
conditions for the matter and metric perturbations in the 
AXCDM model. However, since the ACDM model does 
not include DE perturbations, we should set independent 
initial conditions on Snidi) and Oniai)- As already dis- 
cussed, the scales relevant to the matter power spectrum 
remain always well below the sound horizon ([65]) and we 
expect negligible DE perturbations at any time. Thus, 
the most natural choice for the initial values of the DE 
perturbations is: 

SD{a,)^0, 0D(aO = O. (104) 

Indeed, this is not the only reasonable choice. For in- 
stance, we could have also assumed the adiabatic initial 
condition (f54|) for the DE density contrast, i.e. (^^(ai) = 
{l + Weiai))SM{ai), with ^a/ (a^) given by Eq.(l94|. Again, 
it has been checked that the evolution of the perturba- 
tions does not depend significantly on the particular ini- 
tial condition used. 

Assuming the initial conditions ([M)) . ([55|) and (|104p for 
the matter, metric and DE perturbations at a.^, respec- 
tively, we can solve the perturbed equations (|43|) - (|46)) . 
Equivalently, we can solve . to obtain {Sd^Od) 
and then ((50|) to get the matter density fiuctuations to- 
day SM{k,a = 1) for any dark energy model. In par- 
ticular, we can (as a consistency check) solve the per- 
turbation equations for the ACDM model (in that case 
Soia) — Ooia) — 0, so the only equations needed are ([^5]) 
and (j46[) ). In doing so, we recover exactly the spectrum 
Pa defined in 

Now we proceed to compute the spectrum of the 
AXCDM, PAx{k). In order to better compare the shape 
of the different spectra and the goodness of their fit to 
the 2dFGRS observed galaxy power spectrum, we will 
normalize them at the smallest length scale considered 
£ ~ i.e. at fc = 0.2, taking the ACDM spectrum 

(|96p as reference. To this purpose, we introduce a nor- 
malization factor in the matter power spectrum: 

PAx{k) = AAx\SM{k,a^l)\\ (105) 

Notice that the normalization factor ^ax gives us the 
difference in the matter power spectrum of the model 
with respect to that of the ACDM at the specific scale k — 
0.2. Let us clarify that the reason for choosing this scale 
for the normalization is that, as discussed in section fVl 
the smaller the length scale (i.e. the higher the value of k) 
the less important the DE perturbations are. Therefore, 
at k = 0.2, the matter power spectrum of the model 
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FIG. 3: The ACDM power spectrum Pa(A;) (dot-dashed line) 
versus the normalized and unnormalized spectrum predicted 
by the AXCDM model, for DE sound speeds =0.1 (dashed 
line) and = 1 (solid/gray line): (a) for a set of parameters 
allowed by the analysis of Ref. [23] (of. shaded and striped 
region in our Figs. [321, Q,\ = 0.8, u = uo = 2.6 x 10"^ and 
wx = —1.6. The corresponding curves PA(fc) and PAx{k) 
coincide in this case; (b) for a set of parameters not allowed 
by the F-test fSS] (points in the stripped, but non-shaded, 
region in our Fig. [!}, 0° = +0.35, v = —0.2 and wx = —0.6. 
In this case, PAx(fc) presents a slight deviation as compared 
to PA(fc) at large scales (i.e. at small k). The lower set of 
curves in (b) displays the real (unnormalized) growth, see the 
text. 



should not depend significantly on whether we consider 
the efltect of the DE perturbations; in particular, it should 
be independent of the speed of sound Cg ■ For larger scales, 
however, the DE perturbations can be more significant 
and, as we shall see below, they may introduce some 
differences in the shape of the power spectrum, which 
are nevertheless small in the linear regime. 

The AXCDM power spectrum was calculated for two 
fiducial values of the DE speed of sound, = 1 and = 
0.1 and several combinations of the parameters v, ujx 
and For values of the parameters allowed in Figs.[T][2l 
(shaded and striped region) we find that A\x ~ 1 (within 
~ 10% of accuracy). 

In Fig. [3^, we put together the 2dFGRS observed 
galaxy power spectrum, the ACDM spectrum and the 



normalized and unnormalized AXCDM one, for the set 
of parameters £7^ = 0.8, i/ = i/q = 2.6 x 10~^ and 
Wx = —1.6, which are allowed in Figs.[T][ll For these val- 
ues, we have obtained a normalization factor Aax — 1-1 
and an accurate agreement between the AXCDM power 
spectrum and Pa(^)- This was expected since we are 
assuming allowed values of the parameters, i.e., values 
already consistent with LSS data according to the 'effec- 
tive' approach used in [23| (cf. the discussion in section 
IVII[) . Therefore, the predicted power spectrum from the 
AXCDM ought to be very close to the ACDM one, which 
is in fact what we have substantiated now by explicit nu- 
merical check. 

However, for values of the parameters out of the al- 
lowed region in Fig. [1] the predicted matter power spec- 
trum can differ significantly from the ACDM one, PA(fc). 
This occurs mainly for points that do not satisfy the "F- 
test" condition |23j, even if the other observable con- 
straints (namely the ones related to nucleosynthesis and 
the present value of the EOS (cf. section |VlT|) are ful- 
filled. Let us remind that the F-test consists in requiring 
that the matter power spectrum of the model under con- 
sideration (in this case, the AXCDM model) differs from 
that of the ACDM in less that a 10%, under the assump- 
tion that DE perturbations can be neglected. Given the 
fact (explicitly analyzed here) that the DE perturbations 
should not play a very important role, it is reasonable to 
expect that the F-test should be approximately valid even 
when we do not neglect the DE perturbations. Thus, the 
AXCDM model should exhibit a large deviation in the 
amount of growth with respect to the ACDM precisely 
for those points failing the F-test. Points of this sort 
are those located in the striped region, but outside the 
shaded one in Fig.[T] For these points, we should expect 
an anomalously large normalization factor A\x (namely, 
the factor that controls the matching of the two over- 
all shapes) and, at the same time, we may also observe 
an evident scale dependence in the power spectrum, i.e. 
some significant difference in the predicted shape as com- 
pared to the ACDM one. Such potentially relevant scale 
dependence (or fc-dependence) is introduced by the DE 
perturbations themselves through the last term on the 
r.h.s of ([M)) and is eventually fed into equations 
(gSl). 

As a concrete example, let us consider Fig.[5}D where we 
compare the 2dFGRS observed galaxy power spectrum 
and Pa(A:) with the AXCDM matter power spectrum 
PAx{k) for the following set of parameters: 17° = -1-0.35, 
1/ = —0.2 and wx = —0.6. These values fulfill the nucleo- 
synthesis bound (constraint No. 1 in section lVlH) . specif- 
ically, we have |e| = 0.08 for these parameters (meaning 
that DE density at the nucleosynthesis time represents 
roughly only 8% of the total energy density); and sat- 
isfy also the current EOS constraint (No. 3 in section 
IVII[) : — —0.8. However, this choice of parameters 
largely fails to satisfy the constraint No. 4, i.e. the F- 
test: indeed, we find F = 2.06, which implies that the 
discrepancy in the amount of growth with respect to the 
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ACDM when we neglect DE perturbations is more than 
200%! As expected, for such set of parameters we en- 
counter a large normaUzation factor for the two fiducial 
DE sound speeds = 1 and = 0.1 that we are using 
in our analysis (on average Aax — 2.7). This is reflected 
in the evident gap existing between the upper and lower 
set of curves in Fig. [3]3. The lower set reflects the real 
growth |(5M(fc)P of matter perturbations before applying 
the normalization factor. Such normalization consists in 
the following: for the smallest scale available in the data, 
the AXCDM curves have been shifted upwards until they 
match up with the standard ACDM prediction. Apart 
from the overall gap between the two set of curves, we 
also find a significant shape deviation with respect to the 
standard ACDM model at large scales, as it is patent in 
Fig. This feature is more clearly seen at small sound 
speeds, see next section. 



IX. MATTER AND DARK ENERGY DENSITY 
FLUCTUATIONS 

As we have discussed in section |Vl the DE fluctua- 
tions Sd should oscillate and become eventually negligi- 
ble as compared to the matter fluctuations Sm, specially 
at small scales (inside the sound horizon). However, as 
also noted above, for values that significantly violate the 
F-test 23] , the power spectrum and its shape can be 
noticeably different from that of the ACDM model (cf. 
Fig. [3}d). This suggests that, under such circumstances, 
the DE density perturbations are not completely negligi- 
ble owing to the fact that the term which depends on k 
in the perturbation equations is also proportional to Sd - 
see Eq. (|45|) . In addition, there appears a suppression of 
the growth of matter fluctuations in comparison with the 
growth predicted by the ACDM model. This inhibition of 
matter growth is characteristic of cosmologies where the 
DE behaves quintessence-like, i.e. when the DE density 
decreases with the expansion, whereas phantom-like DE 
(increasing with the expansion) would cause the opposite 
effect (an enhancement of the power). A similar situa- 
tion was also observed in Ref.[4l|] for models with pure 
running A, where in the case v > (in which A decreases 
with the expansion) there is an inhibition of growth while 
for < (when A increases with the expansion) there 
is an enhancement - see also [iO, HH and (g^l for other 
studies. 

Let us clarify that these differences in the amount of 
growth are present even if we neglect the DE perturba- 
tions, see 1235. In fact, the effect of the latter is very 
small, specially for allowed values of the parameters, and 
becomes noticeable only at large scales. At these scales, 
we find that the DE perturbations tend to compensate 
the suppression produced at the background level. This 
slight enhancement is greater the smaller is the DE sound 
speed. Such feature can be appreciated in Fig. 31 where 
we compare the growth of the matter fluctuations at a 
large scale £ ^ (with k = 0.01) predicted by both the 
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FIG. 4: The AXCDM matter density fluctuations at a fixed 
large scale k — 0.01 (in units of /iMpc~^) as a function of 
the scale factor a in comparison with those predicted by the 
ACDM model (dot-dashed line). For the former we have as- 
sumed the same values of the parameters and meaning of the 
lines as in Fig. [S] (a) for the set of parameters in the allowed 
region; (b) for the set of parameters not allowed by the F-test 
in ,23]. 



ACDM and AXCDM models for the two fiducial sound 
speeds of the DE considered before and for the same val- 
ues of the parameters as in Fig. [31 

The growth of matter density fluctuations for the 
AXCDM model is in agreement with the predicted one 
by the ACDM model (the dot-dashed and black line) in 
Fig. IlK, for the set of parameters in the allowed region, 
whereas in Fig.jlla we see the previously commented sup- 
pression for the set of parameters not satisfying the F- 
test. The former case can be compared with Fig. [5^ in 
which we have assumed the same set of allowed parame- 
ters; as expected, we find completely negligible DE fluc- 
tuations today and in the recent past, in agreement with 
the F-test assumption f23| which means completely negli- 
gible DE fluctuations at large scales and a maximum 10% 
of deviation from the ACDM growth of matter density 
fluctuations. On the other hand, values of the parame- 
ters not satisfying the F-test present not only suppression 
on the growth of matter density fluctuations, as shown 
in Fig. Hla, but also larger DE fluctuations today and in 



22 





0.8 t 
1 
1 




0.6 1 






d 


0.4 1 








0.2 [ 




- 




-0? [ 



(a) 



0.2 



0.4 0.6 

a 



0.8 




0.04 0.06 
a 



0.8 
0.6 
0.4 



(b) 



-0.2 



0.2 



0.4 



0.6 



0.8 



FIG. 5: The AXCDM growth of DE fluctuations at the large 
scale k = 0.01 (in units of /iMpc~^) as a function of the scale 
factor a. We have assumed the same values of the parameters 
as in Fig. [3] for DE sound speed =0.1 (dashed line) and 
= 1 (solid/gray line): (a) for the set of parameters in the 
allowed region; (b) for the set of parameters not allowed by 
the F-test in 



the recent past, as shown in Fig. [Sb- 

Furthermore, as discussed in Section |Vl the growth of 
DE fluctuations is expected to oscillate at small scales 
and rapidly decay, what legitimate our assumptions for 
the initial conditions of the DE perturbations. We show 
these oscillations for an allowed set of parameters in 
Fig. [S] Similar behavior is obtained for values of the 
parameters not allowed by the F-test and for both DE 
sound speeds — 1 and = 0.1. The amplitude of the 
DE growth starts negligible 10~^) and rapidly decay 
to zero, as shown in Fig. [S) 

In Fig. [7] we plot the present value of the DE perturba- 
tions as a function of the wave number for the two sets of 
allowed (Fig. [7^) and non-allowed (Fig. [7)d) parameters 
used in the previous plots. We see that the DE pertur- 
bations are negligible at small scales (large k), whereas 
they become larger at larger scales. This is because by 
increasing the scale we are getting closer to the sound 
horizon, as discussed in section |Vl We also see that the 
DE perturbations are larger for the parameters not al- 
lowed by the F-test, which explains why the shape of the 



FIG. 6: The AXCDM growth of DE fluctuations for a small 
scale k — 0.2 (in units of h Mpc~^) and the same set of allowed 
parameters assumed before, and for DE sound speed — 0.1. 



matter power spectrum differs from that of the ACDM in 
this case (cf. Fig. [Sb). However, when comparing with 
Fig. [31 we see that even at the largest explored scale 
(fc = 0.01) the ratio 5d/5m remains rather small, staying 
at the level of 10~^. Finally, let us comment that the DE 
density contrast can become negative with the evolution, 
as in this case happens for = 1. 

As discussed in section |Vl the decay of the DE pertur- 
bations takes place once the term proportional to in 



(I45p becomes dominant. This same term is also responsi- 
ble for the exponential growth of the (DE) perturbations 
in the adiabatic case or, more generally, whenever is 
negative. In order to better appreciate its influence, it is 
useful to compare the evolution of the DE perturbations 
in the adiabatic case (for the most common situation 
where < 0) and the non-adiabatic one (c^ > 0) with 
the scenario in which = 0, since in the latter the term 
proportional to fc^ disappears from the equations. This 
is precisely what has been done in Fig. [5^ for the allowed 
set of parameters used throughout this section. In that 
figure, it is shown the evolution of the DE density con- 
trast at a sufficiently large scale t ^ k^^ for which the 
DE perturbations can be sizeable (namely at A; = 0.01). 
We illustrate the effect for three different regimes of the 
speed of sound: = (with < 0), = and 
= 0.1. The gray lines represent the evolution ol 5d 
when the last term in Ea. (|l^ is neglected, showing in- 
deed that the qualitative behavior of the perturbations 
in the adiabatic and = 0.1 cases does not stem from 
that term. 

For = 0.1, the scale considered is initially (i.e. at 
tti = 1/500) larger than the sound horizon, and thus 
the term proportional to k^ is negligible at the begin- 
ning of the evolution. The same is true for the adia- 
batic case because, in the asymptotic past, the effective 
EOS of the AXCDM model resembles that of matter- 
radiation (u'e(a) Wm for a ^ 0) [l^ and, thus, we 
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FIG. 7: The scale dependence of DE fluctuations today (o = 
1) for the same values of parameters and meaning of the lines 
as in Fig. [5] They rapidly decay at small scales (large k). 
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(We recall that in all our discussion we remain in the 
matter epoch, equality being at a 10^'*). Therefore, 
the term proportional to fc^ is initially unimportant in 
all the three cases and this makes the perturbations to 
evolve in a nearly identical fashion at these first stages, 
as it can be clearly seen from Fig. [5^. 

As the evolution continues, the curve corresponding 
to Cg = 0.1 begin to depart from the others. This oc- 
curs mainly from the instant when the sound horizon is 
crossed, i.e. when the wavelength of the fc-mode gets 
comparable to the sound horizon; such instant can be 
defined through the condition fcA^ = tt, similarly as in 
[4^. Then, the term proportional to fc^ begins to dom- 
inate, which in turn makes the DE velocity gradient 9£i 
to rapidly increase and the DE perturbations to decay. 
Later on, the term proportional to fc^ becomes impor- 
tant also in the adiabatic case. Due to the different sign 
(c^ < 0), the effect that it triggers is now opposite to the 
one observed in the Cg =0.1 case: thus, instead of getting 
stabilized, the DE perturbations initiate an exponential 
growth. 

The previous features can be further assessed in a 
quantitative way by comparing the numerical importance 
(in absolute value) of the two terms inside the curly 
brackets in (M)) . In Fig. [SJs, we plot the ratio between 



FIG. 8: (a) Evolution of the DE perturbations for the same 
set of parameters as in Fig.|3^, at the large scale k = 0.01 
and for three different speeds of sound: = < (solid 
line), = (dashed line) and = 0.1 (dot-dashed line). 
For the latter, the DE perturbations begin to decay after 
the sound horizon crossing (characterized by the condition 
k\s = tt), whereas in the adiabatic case So starts to grow 
exponentially at a ~ 0.2. The last term in Eq. (|44p may be 
neglected (gray lines) without altering the qualitative beha- 
vior, which is triggered by the So-term and, ultimately, by 
the one proportional to k^ in (I45|) : (b) Comparison of the two 
terms inside the curly brackets in (|44p for the =0.1 and 
adiabatic cases. When the So-term becomes important, the 
stabilization (resp. unbounded growth) of the non-adiabatic 
(resp. adiabatic) perturbations becomes manifest. 



these two terms for both the adiabatic case (with < 0) 
and the non-adiabatic situation (c~ = 0.1) (note that the 
term proportional to H^/k^ may be neglected for the 
sub-Hubble perturbations we are dealing with). Com- 
parison with Fig. [5^ reveals that it is precisely when the 
term proportional to 9u stops being negligible that the 
evolution of the perturbations begins to depart from the 
cl = case. The absolute value of the adiabatic speed 
of sound is also shown, in order to illustrate the ultimate 
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reason for the start-up of the exponential growth in the 
adiabatic mode: it is only when begins to depart signif- 
icantly from that the term proportional to k'^ becomes 
important, which in turn triggers a rocket increase of the 
velocity gradient 6o- 

As we have discussed in connection to Fig. [8^, the 
initial evolution of the perturbations in the AXCDM 
model is nearly the same for any of the three values of 
the speed of sound. In fact, in the adiabatic case, and 
given the behavior of the effective EOS in the asymp- 
totic past, the conditions that lead to the simphfied 
setup ([55]) hold. Therefore, that is the equation ini- 
tially controlling the evolution of Sd- For = we 
arrive at exactly the same equation, whereas for positive 

the resulting equation only differs by the last term 
in (j44p . which, as it has been previously discussed, hap- 
pens to be negligible at least during the first stages of 
the evolution. Notice that for We — const., Eq. ([55]) 
integrates to So{a) = (1 -|- We)dM + C, where C is a 
constant determined by the initial conditions. In sec- 
tion UYI we pointed out that C = corresponds to the 
adiabatic initial condition ([5^ . However, for the alter- 
nate initial condition (|104p , and taking into account that 
We{ai) ~ wm = for the AXCDM model in the early 
matter dominated epoch we have C = —5]\i{ai) and 
thus Soia) = 5m {a) ~ 5M{o^i)- From here we find that 
the ratio between DE and matter perturbations in the 
early times of the evolution reads: 



5D{a) 
5 m {a) 



= 1 - 



5M{ai) 
5m (a) 



(106) 



This simple predicted behavior is confirmed from the 
numerical analysis in Fig. [SI where again the allowed 
set of parameters has been used. We see that the ra- 
tio 5d/5m starts being 0, and subsequently as the mat- 
ter perturbations grow the last term in (|106p diminishes, 
until the asymptotic value do /5m = 1 is reached. This 
value is maintained until the conditions leading to (j53p 
cease to be valid. In the adiabatic and > cases, this 
happens when the term proportional to on the r.h.s. of 
Eq. (|45p can no longer be neglected. On the other hand, 
when Cg = 0, the 5d/5m — 1 regime is abandoned at 
the point when the effective EOS starts acquiring sizable 
negative values. Moreover, being the term proportional 
to k^ absent, it is now the last term on the r.h.s. of Eq. 

- which was irrelevant for the other two cases - the 
one that tends to stabilize the DE perturbations (see also 

From the detailed analysis that we have presented here, 
we conclude that the approximation of neglecting the 
DE perturbations can be justified [23[. But this does 
not mean that the computation of these perturbations 
is useless. Indeed, the issue at stake here is not so much 
the quantitative impact of the DE fluctuations upon the 
matter power spectrum - which is actually negligible, as 
we have seen - but rather the fact that the DE perturba- 
tions may be consistently defined in a certain subregion 
of the parameter space only. Of course this subregion 
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FIG. 9: Evolution of the ratio between DE and matter per- 
turbations for the same parameter values as in Fig. |8l The 
special regime of perturbations (|53|l is seen to be approxima- 
tely realized during the first stages of the evolution, for any of 
the three considered speeds of sound. Since We ^ for small 
a in the AXCDM model, we expect So /5m ^ 1 (cf. Eq. (fT06)) ) 
until the conditions leading to (I53|l no longer hold. One can 
clearly confirm this situation in the figure. 



cannot be detected within the context of the effective 
approach [23l |. Therefore, in general, the computation 
of the DE perturbations may have a final quantitative 
bearing on this kind of analyses since it may further re- 
strict the physical region of the parameter space in a very 
significant way. In short, even though the simultaneous 
account of the DE perturbations has a small numerical 
effect on the matter power spectrum within the domain 
where the full system of cosmological perturbations is 
well-defined, it may nevertheless prove to be a highly ef- 
ficient method for excluding large regions of parameter 
space where that system is ill-defined. The upshot is 
that the combined analysis of the DE and matter per- 
turbations may significantly enhance the predictivity of 
the model, as we have indeed illustrated in detail for the 
non-trivial case of the AXCDM model of the cosmic evo- 
lution. 



X. DISCUSSION AND CONCLUSIONS 

In this paper, we have addressed the impact of the 
cosmological perturbations on the coincidence problem. 
In contrast to the previous study [23j , where this pro- 
blem was examined in a simplified "effective approach" 
in which the dark energy (DE) perturbations were ne- 
glected, in the present work we have taken them into ac- 
count in a full-fledged manner. We flnd that the results 
of the previous analysis were reasonable because the DE 
perturbations generally tend to smooth at scales below 
the sound horizon. However, the inclusion of the DE 
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perturbations proved extremely useful to pin down the 
physical region of the parameter space and also to put 
the effective approach within a much larger perspective 
and to set its limitations. 

First of all, we have performed a thorough discussion 
on the coupled set of matter and DE perturbations for 
a general multicomponent fluid. This has prepared the 
ground to treat models in which the DE is a composite 
medium with a variable equation of state (EOS). We have 
concentrated on those cases in which the DE, despite its 
composite nature, is described by a self-conserved den- 
sity p£). Notice that if matter is covariantly conserved, 
the covariant conservation of the DE is mandatory. In 
particular, this is the situation for the standard ACDM 
model, although in this case the self-conservation of the 
DE appears through a trivial cosmological constant term. 
Pa = Pa> which remains imperturbable throughout the 
entire history of the Universe. One may nevertheless en- 
tertain generalized frameworks where the DE is not only 
self-conserved, but is non-trivial and dynamical. This 
is not a mere academic exercise; for instance, in quan- 
tum field theory in curved space-time we generally ex- 
pect that the vacuum energy should be a running quan- 
tity [H, Therefore, in such cases, the CC density 
becomes an effective parameter that may evolve typically 
with the expansion rate, /oa — pa{H), and constitutes a 
part of the full (dynamical) DE of the composite cos- 
mological system with variable EOS. In these circum- 
stances, if the gravitational coupling G is constant, the 
running CC density p\ — p{^{H) cannot be covariantly 
conserved unless other terms in the effective action of 
this system compensate for the CC variation. We have 
called the effective entity that produces such compensa- 
tion "X" or "cosmon" , and denoted with px its energy 
density. Therefore, pD = Pa + Px is the self-conserved 
total DE density in this context, which must be dealt 
with together with the ordinary density of matter pM- 
A generic model of this kind is what we have called the 
AXCDM model [13, Hi. 

Furthermore, from general considerations based on the 
covariance of the effective action of QFT in curved space- 
time m, m. 111] , we expect that the running CC density 
Pa = Ph{H) should be an affine quadratic law of the 
expansion rate i?, see Eq. (IS7)) . Using this guiding prin- 
ciple and the ansatz of self-conservation of the DE, we 
find that the evolution of px , and hence of po , becomes 
completely determined, even though its ultimate nature 
remains unknown. In particular, X is not a scalar field 
in general. 

The AXCDM model was first studied in [l3| as a 
promising solution to the cosmic coincidence problem, 
in the sense that the coincidence ratio r = pd/pm can 
stay relatively constant, meaning that it docs not vary 
in more than one order of magnitude for many Hubble 
times. The main aim of the present paper was to make a 
further step to consolidate such possible solution of the 
coincidence problem, specifically from the analysis of the 
coupled system of matter and DE perturbations. Let us 



remark that this has been a rather non-trivial test for 
the AXCDM model. Indeed, after intersecting the region 
where the DE perturbations of this model can be con- 
sistently defined, with the region where the coincidence 
problem can be solved [TtL |23| | . we end up with a signif- 
icantly more reduced domain of parameter space where 
the model can exist in full compatibility with all known 
cosmological data. The main conclusion of this study is 
that the predictivity of the model has substantially in- 
creased. Therefore, it can be better put to the test in the 
next generation of precision cosmological observations, 
which include the promising DES, SNAP and PLANCK 
projects [63i] . 

Interestingly enough, we have found that the final re- 
gion of the parameter space is a naturalness region which 
is more accessible to the aforementioned precision ex- 
periments. For example, we have obtained the bound 
^ J/ < f ~ 10~^ for the parameter that determines the 
running of the cosmological term. This bound is perfectly 
compatible with the physical interpretation of v from its 
definition Moreover, our analysis indicates that the 
cosmon entity X behaves as "phantom matter" [l7| , i.e. 
it satisfies ujx < — 1 with negative energy density. This 
result is a clear symptom (actually an expected one) from 
its effective nature. It is also a welcome feature; let us re- 
call that "phantom matter" , in contrast to the "stan- 
dard" phantom energy, prevents the Universe from reach- 
ing the Big Rip singularity. Finally, perhaps the most no- 
ticeable (and experimentally accessible) feature that we 
have uncovered from the analysis of the DE perturbations 
in the AXCDM model, is that the overall EOS param- 
eter We associated to the total DE density pu behaves 
effectively as quintessence {we ^ —1) in precisely the 
region of parameter space where the cosmic coincidence 
problem can be solved. In other words, quintessence is 
mimicked by the AXCDM model in that relevant region, 
despite that there is no fundamental quintessence field in 
the present framework. A detailed confrontation of the 
various predictions of the AXCDM model (in particular, 
the kind of dependence We = We{z)) with the future accu- 
rate experimental data [63| , may eventually reveal these 
features and even allow to distinguish this model from 
alternative DE proposals based on fundamental quintes- 
sence fields. 

To summarize: we have demonstrated that the set of 
cosmological models characterized by a composite, and 
covariantly conserved, DE density in which the vac- 
uum energy pA is a dynamical component (specifically, 
one that evolves quadratically with the expansion rate, 
see Eq. (|57| ). proves to be a distinguished class of models 
that may provide a consistent explanation of why pzj is 
near pM, in full compatibility with the theory of cos- 
mological perturbations and the rest of the cosmological 
data. Remarkably, such class of models is suggested by 
the above mentioned renormalization group approach to 
cosmology. We conclude that the AXCDM model can be 
looked upon as a rather predictive framework that may 
offer a robust, and theoretically motivated, dynamical 
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solution to the cosmic coincidence problem. 
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